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ABSTRACT 

MULTIPLE PARTON RADIATION IN 

HADROPRODUCTION 
AT LEPTON-HADRON COLLIDERS 

By 

Pavel M. Nadolsky 

Factorization of long- and short-distance hadronic dynamics in perturbative Quan- 
tum Chromodynamics (QCD) is often obstructed by the coherent partonic radiation, 
which leads to the appearance of large logarithmic terms in coefficients of the pertur- 
bative QCD series. In particular, large remainders from the cancellation of infrared 
singularities distort theoretical predictions for angular distributions of observed prod- 
ucts of hadronic reactions. In several important cases, the predictive power of QCD 
can be restored through summation of large logarithmic terms to all orders of the per- 
turbative expansion. Here I discuss the impact of the the coherent parton radiation 
on semi-inclusive production of hadrons in deep inelastic scattering at lepton-proton 
colliders. Such radiation can be consistently described in the 6-space resummation 
formalism, which was originally developed to improve theoretical description of pro- 
duction of hadrons at e^e" colliders and electroweak vector bosons at hadron-hadron 
colliders. I present the detailed derivation of the resummed cross section and the 
energy flow at the next-to-leading order of perturbative QCD. The theoretical results 
are compared to the experimental data measured at the ep collider HERA. A good 
agreement is found between the theory and experiment in the region of validity of the 



resummation formalism. I argue that semi-inclusive deep inelastic scattering (SIDIS) 
at lepton-hadron colliders offers exceptional opportunities to study coherent parton 
radiation, which are not available yet at colliders of other types. Specifically, SIDIS 
can be used to test the factorization of hard scattering and collinear contributions at 
small values of x and to search for potential crossing symmetry relationships between 
the properties of the coherent radiation in SIDIS, e+e~ hadroproduction and Drell- 
Yan processes. 
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Chapter 1 
Introduction 



Since its foundation in 1970's, perturbative Quantum Chromodynamics (PQCD) has 
evolved into a precise theory of energetic hadronic interactions. The success of the 
QCD theory in the quantitative description of hadronic experimental data originates 
from the following fundamental ideas: 

1. Hadrons are not elementary particles. As it was first shown by the quark model 
of Gell-Mann and Zweig basic properties of the observed low-energy hadronic 
states are explained if hadrons are composed of a few "constituent quarks" with 
spin 1/2, fractional electric charges and new quantum numbers of flavor and 
color 0]. If the hadron constituents {partons) are bound weakly at some energy, 
they can possibly be detected in scattering experiments. The parton model of 
Feynman and Bjorken suggested that the pointlike hadronic constituents may 
reveal themselves in the wide-angle scattering of leptons off hadronic targets [If. 
The first direct experimental proof of the hadronic substructure came from 
the observation of the Bjorken scaling in the electron-proton deep-inelastic 
scattering ||^; subsequently the quantum numbers of partons were tested in a 
variety of experiments 0. 
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2. The elementary partons of QCD are "current" quarks, which interact with one 
another through mediation of non-Abelian gauge fields (gluons) 0. These gauge 
fields are introduced to preserve the local SU{3) symmetry of the quark color 
charges, in accordance with the pioneering work on non-Abelian gauge sym- 
metries by C. N. Yang and R. L. Mills Remarkably, the QCD interactions 
weaken at small distances because of the anti-screening of color charges by 
self-interacting gluons Due to this feature (asymptotic freedom) of QCD, 
probabilities for parton interactions at distance scales smaller than 1 GeV^^ can 
be calculated as a series in the small QCD running coupling as- In the opposite 
limit of large distances, as grows rapidly, so that the QCD interactions be- 
come nonperturbative at the scale of about 0.2 GeV~^. Such scale dependence 
of the QCD coupling explains why the partons behave as quasi-free particles 
when probed in the energetic collision, but eventually are confined in colorless 
hadrons at the later stages of the scattering. 

3. Because of the parton confinement, quantitative calculations within QCD re- 
quire systematic separation of dynamics associated with short and long distance 
scales. The possibility for such separation is proven by factorization theorems 
[p!0| -[T5|. With time, the factorization was proven for observables of increas- 
ing complexity. In 1977, G. Sterman and S. Weinberg introduced a notion of 
infrared-safe observables, which are not sensitive to the details of long-distance 
dynamics |[T^. A typical example of an infrared-safe observable is the cross- 



section for the production of well-separated hadronic jets at an e~^e~ collider. 
It was shown that infrared-safe observables can be systematically described by 
means of PQCD. As a next step, factorization was proven for inclusive observ- 
ables depending on one large momentum scale Q^. In the limit ^ oo, 
such observables can be factorized into a perturbatively calculable hard part. 



describing energetic short-range interactions of hadronic constituents, and sev- 
eral process-independent nonperturbative functions, relevant to the complicated 
strong dynamics at large distances. 

The proof of factorization is more involved for hadronic observables that depend 
on several momentum scales (e.g., differential cross sections). The complications 
stem not so much from the complex dependence of the cross sections on kinematical 
variables, but from the presence of logarithms Inr, where r is some dimensionless 
function of the kinematical parameters of the system. For instance, r may be a ratio 
of two momentum scales Pi and P2 of the system, r = Pi/P2- Near the boundaries 
of the phase space, the ratio r can be very large or very small, in which case the 
convergence of the series in the QCD coupling as can be spoiled by the largeness of 
terms proportional to powers of Inr. Hence the factorization cannot be proven as 
straightforwardly as in the case of the inclusive observables, because its most obvious 
requirement - sufficiently rapid convergence of the perturbative series - is violated. 

To restore the convergence of PQCD, one may have to sum the large logarithmic 
terms through all orders of ag. This procedure is commonly called resummation. 
Logarithmic terms of one rather general class appear due to the QCD radiation along 
the directions of the observed initial- or final-state hadrons (coUinear radiation) or 
the emission of low-energy gluons {soft radiation). Such logarithms commonly affect 
observables sensitive to the angular distribution of the hadrons. In several important 
processes, the soft and collinear logarithms can be consistently resummed through 
the use of the formalism developed by J. Collins, D. Soper and G. Sterman (CSS). 

The original resummation technique was proposed in Ref. fl^ to describe angular 
distributions of back-to-back jets produced at e~^e~ colliders (Fig. |l]l|a) . Subsequent 
developments of this technique and its comparison to the data on the e'^e~ hadropro- 



duction were presented in Refs. 1 18-20 . In Ref. |21 the resummation formalism was 




(b) 

Figure 1.1: (a) Production of hadronic jets at e"^e~colliders; (b) Production of lepton 
pairs at hadron-hadron colliders 
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extended to describe transverse momentum distributions of lepton pairs produced 
at hadron-hadron colliders (Fig.|l.ljb). In the subsequent publications p^-pTfl, this 
technique was developed to a high degree of numerical accuracy. Currently the re- 
summation analysis of this type is employed in the measurements of the mass ||28| and 
the width ||2|| of the PF-bosons produced at the pp collider TEVATRON. With some 
modifications, this resummation formalism is also used to improve PQCD predictions 



for the production of Higgs bosons |^ and photon pairs ||31| at the Large Hadron 
Collider (LHC). 

The hadroproduction at e"'"e~colliders and the lepton pair production at hadron- 
hadron colliders (Drell-Yan process) are the simplest reactions that require resumma- 
tion of the soft and collinear logarithms. In both reactions, the interaction between 
the leptons and two initial- or final-state hadronic systems is mediated by an elec- 
troweak boson V with a timelike momentum. The CSS resummation formalism can 
also be formulated for reactions with the exchange of a spacelike electroweak vector 



boson p2|,|33|. In this work, I discuss resummation in the semi-inclusive production 
of hadrons in electron-hadron deep-inelastic scattering, which is the natural analog of 
e~^e~ hadroproduction and Drell-Yan process in the spacelike channel. The reaction 
of semi-inclusive deep-inelastic scattering (SIDIS) e + A^e + B + X^ where A and 



B are the initial- and final-state hadrons, respectively, is shown in Figure |L2 



As in the other two reactions, in SIDIS the multiple parton radiation affects an- 
gular distributions of the observed hadrons. The study of the resummation in SIDIS 
has several advantages in comparison to the reactions in the timelike channels. First, 
SIDIS is characterized by an obvious asymmetry between the initial and final hadronic 
states, so that the dependence of the multiple parton radiation on the properties of 
the initial state can be distinguished clearly from the dependence on the properties of 
the final state. In contrast, in e"'"e~ hadroproduction or the Drell-Yan process some 
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P B 

Figure 1.2: Semi-inclusive deep inelastic scattering 

details of the dynamics may be hidden due to the symmetry between two external 
hadronic systems. Notably, I will discuss the dependence of the resummed observables 
on the longitudinal variables x and z, which can be tested in SIDIS more directly than 
in e+e" hadroproduction or Drell-Yan process. 

Second, SIDIS can be studied in the kinematical region covered by the measure- 
ments of the hadronic structure functions Fi{x, Q^) in completely inclusive DIS. The 
ongoing DIS experiments at the ep collider HERA probe Fi{x, Q"^) at x down to 10~^, 
which are much smaller than lower values of x reached at the existing hadron-hadron 
colliders. The region of low x, which is currently studied at HERA, will also be 
probed routinely in the production of VF^, and Higgs bosons at the LHC. At such 
low values of x, other dynamical mechanisms may compete with the contributions 
from the soft and collinear radiation described by the CSS formalism. The study of 
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the existing SIDIS data provides a unique opportunity to learn about the applicabil- 
ity of the CSS formalism in the low-x region and estimate robustness of theoretical 
predictions for the electroweak boson production at the LHC. 

Last, but not the least, is the issue of potential symmetry relations between the 
resummed observables in SIDIS, e~^e~ hadroproduction and Drell-Yan process. In 
SIDIS, the dynamics associated with the initial-state radiation may be similar to the 
initial-state dynamics in the Drell-Yan process, while the final-state dynamics may 
be similar to the final-state dynamics in e"'"e~hadroproduction. It is interesting to 
find out if the data support the existence of such crossing symmetry. 



The results presented here were published or accepted for publication in Ref. |35- 



37| . The remainder of the thesis is organized as follows. In Chapter |2|, I discuss 
the basics of factorization of mass singularities in hadronic cross sections. Then I 
review the general properties of the Collins-Soper-Sterman resummation formalism 
and illustrate some of its features with the example of hadroproduction at e"'"e~ 
colliders. 

In Chapter^ I apply the resummation formalism to semi-inclusive deep inelastic 
scattering. Guided by the similarities between SIDIS, e'^e~ hadroproduction and 
Drell-Yan process, I introduce a set of kinematical variables that are particularly 
convenient for the identification and subsequent summation of the soft and collinear 
logarithms. I also identify observables that are directly sensitive to the multiple 
parton radiation. In particular, I argue that such radiation affects the dependence of 
SIDIS cross sections and hadronic energy flow on the polar angle in the photon-proton 
center-of-mass frame. Next I derive the 0{as) cross section and obtain the 0(0:5') 
coefficients for the resummed cross sections and the hadronic energy flow. 

In Chapter^, I compare the results of the CSS resummation formalism and 0{as) 
fixed-order calculation with the data from the ep collider HERA. I show that the CSS 
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resummation improves theoretical description of various aspects of these data. I also 
discuss the dependence of the resummed observables on the longitudinal variables 
X and z. I show that the HERA data are consistent with the rapid increase of 
nonperturbative contributions to the resummed cross section at a; < 0.01. I discuss 
the potential dynamical origin of such low-x behavior of the CSS formula. 

Finally, in Chapter^ I discuss the impact of the multiple parton radiation on az- 
imuthal asymmetries of the SIDIS cross sections. I show that the CSS resummation 
formalism can be used to distinguish reliably between perturbative and nonperturba- 
tive contributions to the azimuthal asymmetries. I also suggest to measure azimuthal 
asymmetries of the transverse energy flow, which provide a clean test of PQCD. 
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Chapter 2 



Overview of the QCD factorization 



Perturbative calculations in Quantum Chromodynamics rely on a systematic proce- 
dure for separation of short- and long-distance dynamics in hadronic observables. The 
proof of feasibility of such procedure naturally leads to the methods for improvement 
of the convergence of the perturbative series when this convergence is degraded by 
infrared singularities of contributing subprocesses. Here I present the basics of the 
factorization procedure. The omitted details can be found in standard textbooks on 



the theory of strong interactions, e.g., Refs. |38-41 



2.1 QCD Lagrangian and renormalization 

Low-energy hadronic states have internal substructure. They are composed of more 
fundamental fermions (quarks) that are bound together by non-Abelian gauge forces. 
The quanta of the QCD gauge fields are called gluons. Quantum ChromoDynamics 
(QCD) is the theory that describes strong interactions between the quarks. In the 
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classical field theory, the QCD Lagrangian density in the coordinate space is 



^qcd{x) = ^ijf{i^-g4aTa-mf)iljf 
f 

- \FfF,^fs-^{r]o,A2y + Ca{SadVd-gCaMVA)cd, (2.1) 
where ipfi{x), A^i^x) and Ca{x) are the quark, gluon and ghost fields, respectively; 

Ff{x) = d^Al - d^Al - gCabcAtAP (2.2) 

is the gauge field tensor; —\{r]a ■ /'^ is the term that fixes the gauge r] ■ A = 
The vector 77° is equal to the gradient vector in covariant gauges {d^A'^ — 0) or an 
arbitrary vector n° in axial gauges {riaA^ — 0). The color indices l,m vary between 
1 to Nc (where A'c = 3 is the number of colors), while the color indices a, b, c, d vary 
between 1 and — 1. The index / denotes the flavor {i.e., the type) of the quarks, 
which is conserved in the strong interactions. The remaining parameters in Cqcd{x) 
are the QCD charge g and the masses of the quarks nif. 

The QCD Lagrangian is invariant under the gauge transformations of the SU {Nc) 
color group: 

M^) ^ u{e{x))ijf{xy, (2.3) 
TaA:{x) ^ u{e{x))TaA:{x)u-\e{x)) + '-{d^u{e{x)))u-\e{x)), (2.4) 

where the x— dependent unitary operator U{9{x)) is 

U{e{x)) = e-^'^»^«W. (2.5) 
{Ta)ijji and Cabd are generator matrices and structure constants of the color group. 
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The commutators of the matrices {Ta)im are 

[Ta, n] = tCabcTc. (2.6) 

The quark fields tpfi and gauge fields A" are vectors in the fundamental and adjoint 
representations of SU (Nc) , respectively. 

In the quantum theory, ipj, A°',Ca are interpreted as "bare" (unrenormalized) op- 
erators of the corresponding fields; g and m/ are interpreted as the "bare" charge 
and masses. The perturbative calculation introduces infinite ultraviolet corrections 
to these quantities. In order to obtain finite theoretical predictions, Cqcd has to be 
expressed in terms of the renormalized parameters, which are related to the "bare" 
parameters through infinite multiplicative renormalizations. 

If the ultraviolet singularities are regularized by the continuation to n = 4 — 2e, 
e > dimensions P2|, the renormalized parameters (marked by the subscript "i?") 



are related to the bare parameters as 

^jnifi) = Z^,\fi),Pj, (2.7) 
AM = Z-,\f,)A:, (2.8) 

Cai?(/i) = Z^^{fl)Ca, (2.9) 

gM = Zg\f,)fi-^g, (2.10) 
mfnifi) = Z;^^{ii)mj, (2.11) 

where Z^, Za, Z^^ Zg, and Zm are perturbatively calculable renormalization constants. 
In the dimensional regularization, the renormalized parameters depend on an auxil- 
iary momentum scale /i„, which is introduced to keep the charge g dimensionless in 
7^ 4 dimensions. In Eqs. p.7| - |2.11| ) the renormalized parameters and the constants 
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are expressed in terms of another scale which is related to //„ as 



(2.12) 



Here 7^ = 0.577215... is the Euler constant. 

2.2 Asymptotic freedom 

The further improvement of the theory predictions for physical observables is achieved 
by enforcing their invariance under variations of the scale i.e., by solving renormal- 
ization group (RG) equations. Consider an observable S that depends on N external 
momenta pf , i — 1, . . . , N. If the renormalized expression for S is 

S {9R{fj), {m/ij(/x)}, {pi}, /i) 

(where "{...}" denotes a set of parameters), then the RG-improved expression for S 
is 



where g{iJi) and rufdj,) are the running QCD charge and quark masses. By solving 
the equation for the independence of ^S" from /x, 



(2.13) 



(2.14) 



we find the following differential equations for g{ii) and m/(//): 
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= m^^)), (2.15) 

^^a^ = -lmf{g{fi))mf{fi). (2.16) 



The approximate expressions for the functions P{g) and ■jm{g) on the r.h.s. of 
Eqs. (|2.15| ) and (|2.16|) are found from the /x— dependence of the fixed-order renor- 



malized charges and masses: 

Pignif.)) ^ /.^, (2.17) 
1 dm%{fi) 

The renormalization group analysis of the QCD Lagrangian suggests that the 
interactions between the quarks weaken at high energies, i.e., that Quantum Chro- 
modynamics is asymptotically free in this limit. Indeed, the perturbative series for 
the function P{g) is 



oo , 



(^(9) = -gJ2[^) (2-19) 

k=i 



where as = 5'^/47r is the QCD coupling. In the modified minimal subtraction (MS) 
regularization scheme ||4^, the lowest-order coefficient Pi in Eq. (|2.19|) is given by 



Pi = Y^A-^TuNf, (2.20) 
where Nf is the number of active quark fiavors, Ca = Nc = 3, and Tr = 1/2. By 
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solving Eq. (|2.15|) , we find that 

(^s{^A = . . ^ , „2 ■ (2.21) 



This equation proves the asymptotic freedom of QCD interactions: for six known 
quark generations, /3i > and 



hm asifJ') = 0. 



Higher-order corrections to the beta-function do not change this asymptotic behavior. 
Eq. ( |2.21|) also shows that Q;s(/i) has a pole at some small value of /i. The position of 
this pole can be easily found from the alternative form of Eq. (|2.21| ), 

Air 



In Eq. ( 2.22|) , Aqqd is a phenomenological parameter, which is found from the analysis 



of the experimental data. The most recent world average value of Aqcd for Nf = 5 
and expression for the /5-function is 208+1 MeV ||4||. According to Eq. (|2]2|), 



as{^) becomes infinite when /i = Aqcd- This feature of the QCD running coupling 
obstructs theoretical calculations for hadronic interactions at low energies. 



2.3 Infrared safety 

Due to the asymptotic freedom, the calculation of QCD observables at large /i can 
be organized as a series in powers of the small parameter g{n). To find out when the 
perturbative calculation may converge rapidly, consider the formal expansion of the 
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RG-improved expression (|2.13|) for the observable 5* in the series of 



S = *(te}, (2.23) 

In this expression, the function {mj},/i) includes all coefficients that do not 

depend on the order of the perturbative calculation (for instance, the phase space 
factors). The mass dimension of ^{{pi}, {m/},/i) is equal to the mass dimension of 
S. The sum over k on the right-hand side is dimensionless. The coefficients of the 
perturbative expansion S^'^^ depend on dimensionless Lorentz-invariant combinations 
of the external momenta pf, the mass parameter /i, and the running quark masses 



rfif{ij). There are indications that the perturbative series in Eq. (|2.23| ) are asymptotic 
45|| , so that it diverges at sufficiently large k. However, the lowest few terms of this 
series may approximate 5* sufficiently well if they do not grow rapidly when k increases. 
The factors that control the convergence of Eq. (|2.23|) can be understood in a 



simpler case, when all Lorentz scalars Pi ■ pj in Eq. (|2.23| ) are of the same order Q^. 
Then Eq. (|2]2|) simplifies to 



S = $(te}, {mj}, /.) g S(^^ {^^}) fM- (2.24) 

When ^ ^qcdi can choose /i ~ Q to make g{fi) small. This choice also elim- 
inates potentially large terms like ln((5^//U^) from the coefficients S''''\ In addition, 
let's assume that Q is much larger than any quark mass frif{n) on which S depends. 
For instance, 5* may be dominated by contributions from the m, d, s quarks, whose 
running masses are lighter than 200 MeV at /i = 2 GeV At /i' > 2 GeV, the 
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quarks become even lighter due to the running oimj: 

rhfifi') = Thfifi) expl- J y7m/(^) > < %(/^), (2.25) 



since in QCD 



7n^/(/i) = ^^^Cp + 0{al) > 0. (2.26) 



Here Cp = {N^ - 1)/{2N^) = 4/3. 

When the quark masses vanish, many observables, which are finite if m^j ^ 0, 
acquire infrared singularities. These singularities are generated from the terms in 
the perturbative coefficients that are proportional to the logarithms ln(mj//i^). The 
expansion in the perturbative series ( |2.24|) makes sense only for those observables S 



that remain finite when mflfiY / fi^ — > 0. 

There are two categories of observables for which the perturbative expansion ( p. 241 ) 



is useful. In the first case, the coefficients S^'^^ are finite and analytically calculable 
when /i — > oo : 

«ijm)--(S)-({(^)"}). «>M-) 

Such observables are called infrared-safe ||1^. For instance, the total and jet pro- 
duction cross sections in e~^e~ hadroproduction are infrared-safe. In this exam- 
ple, hadrons appear only in the completely inclusive final state. According to the 
Kinoshita-Lee-Nauenberg (KLN) theorem [H], such inclusive states are free of in- 



frared singularities, so that the finite expressions for the total and jet cross sections 
can be found from the massless perturbative calculation. 
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Figure 2.1: Factorization of collinear singularities in completely inclusive electron- 
hadron DIS 

In the second case, S^^^ are not infrared-safe, but all mass singularities of S^^^ can 
be absorbed (factorized) into one or several process-independent functions. These 
functions can be measured in one set of experiments and then used to make predictions 
for other experiments. 

To understand which singularities should be factorized, notice that there are two 
classes of the infrared singularities in a massless gauge theory: soft singularities and 
collinear singularities. The soft singularities occur in individual Feynman diagrams 
when the momentum k'^ carried by some gluon line vanishes {k'^ ~ Xk'^, where A — > 
and are fixed). The soft singularities cancel at each order of 0:5 (/x) once all Feynman 
diagrams of this order are summed over. 

In contrast, the collinear singularities occur when the momenta pi and P2 of two 
massless particles are collinear to one another , i.e., when pi • P2 — ^ 0. Since one or 
both collinear particles can be simultaneously soft, the class of the collinear singu- 
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larities partially overlaps with the class of the soft singularities. The soft collinear 
singularities cancel in the complete fixed-order result just as all soft singularities 
do. On the contrary, the singularities due to the collinearity of the particles with 
non-vanishing momenta do not cancel and should be absorbed in the long-distance 
phenomenological functions. 

As an illustration of the factorization of the purely collinear singularities, consider 
the factorized form for the cross section of inclusive deep inelastic scattering e + A — > 
e + X (where A is a hadron) in the limit — > oo: 

This representation and notations for the particle momenta are illustrated in Fig. pTT| . 

In Eq. ( p.28| ), = —q^ is the large invariant mass of the virtual photon 7*, 
X = Q"^ /{2{pA ■ Q'))- These variables are discussed in more detail in Subsection p. 
^hard j rfQ^) is the iufrarcd-safc ("hard") part of the cross-section for the scattering 
e + a— >e + Xof the electron on a parton a. Fa/A{iai IJ^'f) is the parton distribution 
function (PDF), which absorbs the collinear singularities subtracted from the full 
parton-level cross section to obtain rfa^"^'^/ {dx dQ"^). In the inclusive DIS, all collinear 
singularities appear due to the radiation of massless partons along the direction of 
the initial-state hadron A. The final state is completely inclusive; hence, by the KLN 
theorem, it is finite. 

The collinear radiation in the initial state depends only on the types of a and A 
and does not depend on the type of the particle reaction. Therefore, Fa/A{^a, fJ^p) 
is process-independent. It can be interpreted as a probability of finding a massless 
parton a with the momentum C,aPA the initial hadron with the momentum p^. To 
obtain the complete hadron-level cross section, we sum over all possible types of a 
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(a = g,u,u,d,d, . . .) and integrate over the allowed range of the momentum fractions 

In Eq. (|2.28|) , both the "hard" cross sections da'^"-'"^ / {dx dQ'^) and the parton dis- 
tribution functions Fa/A{^a, I^f) depend on an arbitrary factorization scale /ij?, which 
appears due to some freedom in the separation of the collinear contributions included 
in Fa/A{^a, ^J'f) from the "hard" contributions included in da^"-^'^ / {dx dQ^) . Of course, 
the complete hadron-level cross section on the l.h.s. of Eq. ( |2.28D should not depend 
on hf- Hence the /Xir-dependence of Fa/ A{^a, f^'p) should cancel the /ii?-dependence of 
the hard cross section. This requirement is used to find Dokshitser-Gribov-Lipatov- 



Altarelli-Parisi (DGLAP) differential equations which describe the dependence 

of Fa/A{^a, fJ'P) on fip- 



fiF- 



dFa/Ai^a, /^f) 



dni 



$^(^f5®^6M)(ea,/iF). (2.29) 



Here are "spacelike" splitting functions that are currently known up to 

0{ag) [Q. They describe evolution of partons with spacelike momenta. The convo- 
lution in Eq. (|2.29| ) is defined as 

{f^g){x,fi)= [ f{x/^,fi)g{^,fi)^. (2.30) 



A similar approach can be used to derive factorized cross sections for reactions with 
observed outgoing hadrons. Such cross sections depend on fragmentation functions 
(FFs) Ds/bi^, f^o), which absorb the singularities due to the collinear radiation in 
the final state. The fragmentation function can be interpreted as the probability of 
finding the hadron B among the products of fragmentation of the parton b. The 
variable ^ is the fraction of the momentum of b that is carried hy B. In the presence 
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of FFs, the hadron-level cross section becomes dependent on yet another factorization 
scale jiD- Similarly to the PDFs, the dependence of the FFs on hd is described by 
the DGLAP evolution equations: 



where are the "timelike" splitting functions. 

As in the case of the renormalization scale /x, it is natural to choose the factor- 
ization scales and fin of order Q to avoid the appearance of the potentially large 
logarithms \n{Q/^p) and \n{Q/fi£)) in the "hard" cross section. I should emphasize 
that the factorized cross sections are derived under the assumption that all Lorentz 
scalars pi ■ pj are of order Q^, so that x in Eqs. ( p.28| ) is sufficiently close to unity. 
When some scalar product Pi ■ pj is much larger or smaller than Q^, the convergence 
of the perturbative series for the hard cross section is worsened due to the large loga- 
rithms of the ratio Pi -Pj/Q"^. This is a general observation that applies to any PQCD 
calculation. In some cases, the predictive power of the theory can be restored by the 
summation of the large logarithms through all orders of the perturbative expansion. 
In particular, the resummation of the large logarithms is required for the accurate 
description of the angular distributions of the final-state particles, including angular 
distributions of the final-state hadrons in SIDIS. In the next Section, I discuss general 
features of such resummation on the example of angular distributions of the jets in 
e'^e~ hadroproduction. 




(2.31) 
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Production of hadrons in decays 
a) No QCD radiation 




e 



q 

b) QCD radiation 




Figure 2.2: The space-time picture of hadroproduction at e~^e colliders 
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2.4 Two-scale problems 



2.4.1 Resummation of soft and collinear logarithms 

To understand the nature of the problem, consider the process e~^e~ ^ ^ jets 
(Fig.p..l|a). The space-time picture of this process is shown in Figure Let us 
assume that the Z°-bosons are produced at the resonance {Ee+ + E^- = Mz) at rest 
in the laboratory frame. In e"'"e~ hadroproduction, the hadronic decays are initiated 
predominantly by the direct decay of the Z^-boson into a quark-antiquark pair. The 
QCD radiation off the quarks produces hadronic jets, which are registered in the 
detector. 

If no additional hard QCD radiation is present (Fig.|2.2|a), the decay of the 
boson produces two narrow jets escaping in the opposite directions in the lab frame. 
The typical angular width of each jet is of the order Kqcd/ Ea,b 1, where Ea,b ~ 
Mzjl are the energies of the jets. The quarks may also emit energetic gluons, in 
which case the angle between the jets is not equal to tt (Fig. p.2|b). If the angle d 



in Fig. p.2|b is large, the additional QCD radiation is described well by the rapidly 
converging series in the small perturbative parameter[| as{Mz)lT^ ■ But when — 0, 
the higher-order radiation is no longer suppressed, because the smallness of as(M^)/7r 
is compensated by large terms ^^(6*^/4)/^^, p > in the hard part of the hadronic 
cross section. Therefore, the calculation at any fixed order does not describe reliably 
the shape of the hadronic cross section when 6' — > 0. 



To illustrate this point, consider the hadronic energy- energy correlation [|49| , de- 
fined as 

dEA / dEB EaEb—pt-tf^ n- (2-32) 



dcos9 My Jq Jq dEAdEsd cos 9 



*From here on, I drop the "bar" in the notation of the running as- 
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In the limit 6* — > 0, but 6' 7^ , dT,/d cos 9 behaves as 

k 2k-l ,„2~ 







dcos9 





fc=l ^ ^ m=0 ^ ^ 



where c^m are calculable dimensionless coefficients. Additionally there are virtual 
corrections to the lowest order cross section, which contribute at = 0. Suppose we 
truncate the perturbative series in Eq. ( p.33| ) at k = N. UN increases by 1 (that 



is, if we go to one higher order in the series of as), the highest possible power of 
the logarithms ln'"(6'^/4) on the r.h.s. of Eq. ( p.33| ) increases by 2. Therefore, the 
theoretical prediction does not become more accurate if the order of the perturbative 
calculation increases. Equivalently, the energy-energy correlation receives sizeable 
contributions from arbitrarily high orders of as- 

To expose the two-scale nature of this problem, let us introduce a spacelike four- 
vector and a momentum scale qx as 



Qt = Q^-Pa Pb > 2.34 

Pa - Pb Pa- Pb 

g2 = -g,%^>0, (2.35) 



where q^, p^, are the momenta of the Z°-boson and two jets, respectively. The 
vector is interpreted as the component of the four-momentum q^ of the Z'^-boson 
that is transverse to the four-momenta of the jets; that is, 

qfPA = qt-PB = 0. (2.36) 

The orthogonality of to both and follows immediately from its defini- 
tion (PI). 
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In the laboratory frame, 



= (M^, 0); 
Pa = Ea{1, ua); 
Pb = Eb{1, -ub), 



(2.37) 
(2.38) 
(2.39) 



where EajUa and Eb, —ub are the energies of the jets and the unity vectors in the 
directions of the jets, respectively. The large invariant mass — M| of the Z°-boson 
can be associated with the QCD renormalization scale Q^. Let the z-axis be directed 
along Ha- Then qT coincides with the length of the transverse component (Jt of : 



9 6 
q^ = ( -Mz tan -, qr, 0, -Mz tan - ) . 



At the same time 



1 — cos 9 



g2 1 + COS^' 



(2.40) 



and 



hm -^4t = — IH h . . . 

e-oQ2 4 V 6 



(2.41) 



We see that the problems at ^ — > arise due to the large logarithmic terms 

In"^(g2/g2)/g2 ^^^^ ^2/g2 ^ 1. 



dcosO 



g2 dE 



2 d^f. 



«5(Q) 



TT 



2/c-l 



m=0 



m , St 



(2.42) 
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Figure 2.3: The structure of infrared singularities in a cut diagram D for the energy- 
energy correlation in the axial gauge 



The origin of these logarithms can be traced back to the presence of infrared 
singularities in the QCD theory. Before considering these singularities, notice that 
the energy-energy correlation is sufficiently inclusive to be infrared-safe. Therefore, 
the complete expression for the energy-energy correlation is finite at each order of 
0:5 (/i). On the other hand, the infrared singularities do appear in individual Feynman 
diagrams. According to the discussion in Section |2]^, these singularities are due to 
the emission of soft gluons.p| Although the soft singularities cancel in the sum of all 
Feynman diagrams at the given order of as, this cancellation leaves large remainders 



where 



(2.43) 



■km ~ A (^km- 



W{ql/Q^)/ql ifqr is small. 



^The purely collinear singularities do not appear because of the overall infrared safety of the 
energy-energy correlation. 



25 



Fortunately, not all coefficients cj,^ in Eq. (|2.42|) are independent. Refs. ||50|,|51 



suggested that the leading logarithmic subseries in Eq. ( |2.42|) and in analogous ex- 
pressions in SIDIS and Drell-Yan process can be summed through all orders of 0:5. 



The possibility to sum all logarithmic subseries in Eq. ( |2.42|) and restore the conver- 
gence of the series in as was proven by J. Collins and D. Soper Schematically, 
Eq. ( ^.42|) can be written as 



as{L + l) 
+ al{L^ + + L + 1) 
+ al {1^ + 1^ + ^ + 1^ + 1 + I) 

. ...}, 



(2.44) 



where L = \n{qj./Q ), and the coefficients 2c'^^/(7r Q ) are not shown. This series 
can be reorganized as 



dq? 



T 



-J {"5^1 + tt|^2 + • • • } , 



where 



asZi ~ asiL + 1) + aliL^ + L^) + aliL^ + L^) + . . . 

aUL + l) + aUL^ + L^) + ... 

al{L + l) + ... 



2 ry 

asZ2 





Bu 


Co 


A2, 


B2, 


Ci 


As 


B3 


C2 



(2.45) 



In Eq. ( |2.45|) , the right-hand side shows the new coefficients AkyB^Ck-i that are 
required to calculate each new subseries agZ^. The complete subseries a^Zk can be 
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reconstructed as soon as the coefficients A^, Bk Ck-i are known from the calculation 
of the term Oi%{L + 1). Each successive subseries a'gZk in Eq. ( |2.45| ) is smaller by 



as than its predecessor, so that as regains its role of the small parameter of the 
perturbative expansion. 

The rule that makes the resummation of the subseries possible follows from 

(a) the analysis of the structure of the infrared singularities in the contributing Feyn- 
man diagrams at any order of as{lJi) and (b) the requirement that the full energy- 
energy correlation is infrared-safe and gauge- and renormalization-group invariant. 

The structure of the infrared singularities can be identified from the analysis of 
analytic properties of the Feynman diagrams with the help of the Landau equations 
||5^H54| and the infrared power counting ^]. This structure for some con- 



tributing cut diagram D is illustrated by Figure p73|. Throughout this discussion the 
axial gauge C ■ ^ = is used.0 In D we can identify two jet parts Ja^ Jb, the hard 
vertex H, and possibly the soft subdiagram Sg. By their definition, the jet parts J a 
or Jb are the connected subdiagrams of D that describe the propagation of nearly 
on-shell massless particles inside the observed jets. Each of the particles in the jet 
part J A has a four-momentum that is proportional to the momentum of the jet A: 

Pt = P^P% (2.46) 
where 

0<A<1, and ^[3i = l. (2.47) 

i 

Similar relations hold for the momenta of the particles in the jet part Jb- 

^The discussion of the infrared singularities in covariant gauges can be found, for instance, in 
Ref.rtsl. 
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Figure 2.4: Examples of the finite soft subdiagrams: (a) the subdiagrams that are 
connected to Ja, Jb by one or several quark lines; (b) the subdiagrams that are 
connected to H 



Both jets originate from the hard vertex H that contains contributions from the 
highly off-shell particles. In the axial gauge the jet parts are connected to H only 
through the single quark lines. Since the hard scattering happens practically at one 
point, H depends only on and not on q^. 

After the jets are created, they propagate in different directions with the speed 
of light. Due to the Heisenberg uncertainty principle, these jets, which are separated 
by large distances, do not interact with one another except by the exchange of low 
momentum (soft) particles. The infrared singularities, which are associated with the 
long-distance dynamics, can occur only in the jet parts or the soft subdiagram. This 
observation can be refined by the dimensional analysis of the Feynman integrals in 
the infrared limit [ij, |3^, |5^ , which shows that the infrared singularities are at most 
logarithmic. Also, those soft subdiagrams that are attached to the jet parts J a and 
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Jb with one or more quark propagators (Fig.|2]^a) or are connected to H (Fig.|2?^) 
are finite. 

To summarize, the infrared singularities of any individual Feynman diagram reside 
in the soft parts of the jets J4, Jb and in subdiagrams Sg that are connected to J4, Jb 
by soft gluon lines (cf. Fig.|2.3|). Both types of singularities contribute at = (i.e., 

= p^j^ +Pb)i ill agreement with the expectation that the small-gr logarithms are 
remainders from the cancellation of such singularities. Therefore, at small qt the 
distribution dL/dq^ naturally factorizes as 



where H{Q'^) is the contribution from the pointlike hard part, S{q'^,Q^) is the all- 
order sum of the large logarithms, and C^^^ collect finite contributions from the jet 
parts. Clearly, C^"* = C^"* due to the symmetry between the jets. 

The factorized formula is proven by considering the Fourier-Bessel transform of 
dT,/dq^ to the space of the impact parameter b, which is conjugate to qx- Explicitly, 



i7(g2)crcr5(g|,g2). 



(2.48) 




(2.49) 



where 



W^s(&,Q) = $^?5cr*(Ci,C2)C^"*(Ci,C2)e-^(^'«'^^'^^). 



(2.50) 



j 



In Eqs. (|2.49| , |2.5U| ) 



(2.51) 
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is the square of the center-of-mass energy of the initial-state electron and positron; 

denotes the summation over the active quark flavors {i.e., j = u,u,d,d, . . .); 
Cj are the couplings of the quarks to the Z^-bosonsI; (o'o/'S'e+e-) ^ the Born 
approximation for the hard part Hi^Q"^). The Fourier-Bessel transform of the shape 
factor S{q^,Q^) is given by e-^C'-^'^i-'^^)^ where S{b,Q,Ci,C2) is called the Sudakov 
function. At 6^ <^ ^qcd {i-^-, in the region of applicability of perturbative QCD), 
the Sudakov function is given by the integral between two momentum scales of the 
order Q and 1/6, respectively: 

limS(6,g,Ci,C2) = / ^L4(a5(7l),Ci)ln^ + S(a5(7l),Ci,C2) , 

(2.52) 

where A and B can be calculated in PQCD. C\ and C2 are arbitrary constants of the 
order 1 that determine the range of the integration in S{b,Q). The undetermined 
values of these constants reflect certain freedom in separation of the collinear-soft 
contributions included in S^{b, Q) from the purely collinear contributions included in 
C-functions. At each order of as, changes in S^{b, Q) due to the variation of Ci, C2 
are compensated by the opposite changes in the C-functions. Hence the perturba- 
tive expansion of W^, does not depend on these constants. However, the complete 
form- factor VFs in Eq. ( p. 501) does have residual dependence on Ci, C2 because of the 
exponentiation of the terms depending on Ci and C2 in exp (— S'(6, Q, Ci, C2)). The 

variation of Ci, C2 allows us to test the scale invariance of the separation of soft and 
§For the up quarks, 

where e is the charge of the positron and 9w is the weak mixing angle. For the down quarks, 



30 



collinear contributions in the Wj^-tevm. 

At 6^ ^ ^QCD^ the behavior of S is determined by complicated nonperturbative 
dynamics, which remains intractable at the current level of the development of the 
theory. At large b the Sudakov function S is parametrized by a phenomenological 
function S^^{b, Q), which has to be found from the comparison with the experimental 
data. When Q ^ oo, the sensitivity of the resummation formula to the nonpertur- 
bative part of S{b, Q) is expected to decrease. 

Now suppose that the experiment identifies a hadron Ha in the jet J a and a hadron 
B in the jet Jb- Let za,b be the fractions of the energies of the jets J a and Js carried 
by Ha and Hb, respectively. The cross section of the process e+e^ — > HaHbX is 
no longer infrared-safe because of the collinear singularities due to the fragmentation 
into the hadrons Ha and Hb. Nonetheless, in the limit qx ^ the cross section 
da / {dzAdzsdq^) factorizes similarly to Eqs. (|2.49|J2.50| ): 



da 



HaHb 



dzAdzBdq, 
where at 6 — ^ 



" J^2'^"'"''W»-«J'>-^-*,zb), (2.53) 



Wh^Hb (b, za, zb) = ^ej 

j 

a,b = g, u , d 

3 =u\ ^d\ .... (2.54) 

The only major difference between the form-factor WhaHb the hadron pair pro- 
duction cross section da / {dzAdzsdq^) and the form-factor H^s for the energy-energy 



31 



correlation dL/dq^ is the presence of the fragmentation functions Dn/aiCy f^), which 
absorb the collinear singularities due to the final-state fragmentation into the observed 
hadrons Ha,Hb. The FFs are convolved with the coefficient functions 
^ab^i^^ Ci,C2, ^J'D,b), which absorb finite contributions due to the perturbative collinear 
radiation. 

The same resummation technique can also be applied to the production of vector 
bosons {e.g., virtual photons 7*, which decay into lepton-antilepton pairs) at hadron- 
hadron colliders (Fig.|l.l|b). In this process, the four- vector is introduced using 
the same definition (|2.34| ), where now and pb denote the momenta of the initial 
hadrons A and B. The scale qr is just the magnitude of the transverse momentum pr 
of 7* in the center-of-mass frame of the hadron beams (Fig. |2.5|) , since in this frame 

= (0,pT,0,0). 

Therefore, the 6-space resummation formalism ||2l| applies to the production of vector 
bosons with small transverse momenta. The cross section for the production of the 
virtual photon 7* at can be factorized as 



dcFj 



dQ'^dydq^ 



T 



/• d% (2.55) 



g^^O ^AB J (27r 



where and y are the virtuality and rapidity of 7* in the lab frame, xa,b = 
and 



(2.56) 



6^0 . 



In W^, Cj are fractional electric charges of the quarks (e^ = 2/3 for up quarks and 
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—1/3 for the down quarks). (Cj" ® Fa/A) {xa, b, jj) and (Cj^ ® -^fe/s) (a^s, At) are the 
jet parts corresponding to the incoming hadrons A and B. They are constructed from 
the perturbatively calculable coefficient functions C^abi^^ f^) convolved with the PDFs 
for the relevant partons. The perturbative part of the Sudakov function in Eq. ( p.56| ) 
has the same functional dependence as in Eq. (|2.52| ) for e"'"e~-hadroproduction. As in 
the case of VFs and WhaHb i the large-6 behavior of should be parametrized by a 
phenomenological function. 

To conclude, the &— space resummation formalism was originally derived to de- 
scribe the production of hadrons at e+e~ colliders fl^ and production of electroweak 



vector bosons at hadron-hadron colliders |^. The possibility to apply the same 
formalism to SIDIS relies on close similarities between the three processes. First, 
hadronic interactions in all three processes are described by the same set of Feynman 
diagrams in different crossing channels. Second, multiple parton radiation dominates 
each of the three processes when the final-state particle escapes closely to the direc- 
tion predicted by the leading-order kinematics. The formalism for the resummation 
of such radiation can be formulated in Lorentz-invariant notations, so that it can be 
continued from one process to another. 
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Vector boson production at 
hadron-hadron colliders 




b) QCD radiation 

Pt ^ 




Figure 2.5: The space-time picture of Drell-Yan process 



34 




Figure 2.6: The ladder structure of the DIS cut diagrams 



2.4.2 QCD at small x 



According to the discussion in Section |2]^, the convergence of the series in as{ij) 



depends on the absence of very large or small dimensionless quantities in the pertur- 
bative coefficients. In particular, the dimensionless variable x in the inclusive DIS 
cross section should not be too close to zero: otherwise the hard part of the DIS cross 
section contains large logarithms ln'^(l/x), which compensate for the smallness of 
as{Q)- These logarithms are different from the logarithms ln™/i^ resummed by the 
DGLAP evolution equations. As a result, the factorization of the DIS cross section in 
the hard cross section and PDFs (cf. Eq. (|2.28|) ) may experience difficulties at small x. 

The large logarithms ln'"(l/x) are resummed in the formalism of Balitsky, Fadin, 
Kuraev and Lipatov (BFKL) ||5^. The BFKL and DGLAP pictures for the history of 
the parton probed in the hard scattering are quite different. Both types of formalisms 
resum contributions from the cut ladder diagrams shown in Figure |2^ . In this Fig- 
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ure, each "rung" Rk is a two-particle irreducible subdiagram that corresponds to the 



radiation off the probed parton line (see Refs. ||3^,|53l for more details). The vertical 
propagators correspond to the quarks or the gluons that are parents to the probed 
quark. The momenta 2 n ^'^^ from the parent hadron to the probed quark. The 
momenta 2,...,n Ao'^ through the rungs and are sums of the momenta of the radiated 
particles. The conservation of the momentum in each rung implies that 

Pi = pf+i - fcf , i = 1, . . . , n, (2.57) 

where Pn+i = Pa- the reference frame where the hadron A moves at the speed of 
light along the z-axis, x coincides with the ratio of the plus components of ^ind p^: 



x = ^, (2.58) 
Pa 



where 



k^ = k°± k\ (2.59) 

The DGLAP equation arises from the resummation of the ladder diagrams corre- 
sponding to the collinear radiation along the direction of the hadron A. The radiating 
parton remains highly boosted at each rung of the ladder. At the same time, the trans- 
verse momentum carried away by the radiation grows rapidly from the bottom to the 
top of the ladder. The DGLAP equation corresponds to the strong ordering of the 
transverse momenta flowing through the rungs Ri'. that is, 

> 4i > ^2 > • • ■ > 4„ > Aqcd, (2-60) 
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while 



Pl-Pt Pn --Pa>0 (2.61) 



and 



Pi ~ P2 Pn ~ Pa ~ 0- (2-62) 

On the other hand, the BFKL formalism describes the situation in which the QCD 
radiation carries away practically all energy of the probed parton. In this case, 

< < • • • < p+ < p1, (2.63) 

and 

Pr > P2 > • • • > Pn > p:4- (2-64) 

In addition, the BFKL picture imposes no ordering on the transverse components 
of k^: 



4l ^2 4n » ^QCD- (2.65) 



As a result, the probed quark is likely to have a significant transverse momentum 
throughout the whole process of evolution, which is impossible in the DGLAP picture. 
Due to its large kx, the radiating parton is olf its mass shell at any moment of its 
evolution history, so that the BFKL radiation cannot be factorized from the hard 
scattering. As another consequence of the ^T-unordered radiation, the BFKL picture 
implies broad angular distributions of the final-state hadrons, while in the DGLAP 
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picture the hadrons are more likely to belong to the initial- and final-state jets. 

Since the BFKL approach applies to the limit x — and kfj^ ^ ^qcdi it corre- 
sponds to asymptotically high energies of hadronic collisions. So far, the experiments 
have produced no data that would definitely require the BFKL formalism to explain 
them. In particular, the behavior of the inclusive DIS structure functions in the 
low X region at HERA agrees well with the predictions of the traditional fac- 

torized formalism and disagrees with the steep power-law growth predicted by the 



leading-order solution of the BFKL equation ||56 



The situation is not so clear for some less inclusive observables, which deviate 
from the low-order predictions of PQCD. Specifically, SIDIS in the small-x region 
is characterized by large higher-order corrections. Some of these corrections can be 
potentially attributed to the enhanced /cr-unordered radiation at x — 0. If this is 
indeed the case, the effects of the fcr-unordered radiation may be identified by observ- 
ing the changes in the angular distributions of the final-state hadrons or "intrinsic ^t" 
of the partons. In order to pinpoint these effects, good understanding of the angular 
dynamics in the traditional DGLAP picture is needed. Such understanding can be 
achieved in the framework of the small-g-r resummation formalism, which systemati- 
cally describes angular distributions of the hard, soft and collinear radiation. Hence 
it can be naturally used to organize our knowledge about the angular patterns of the 
DGLAP radiation and search for the effects from new \ow-x QCD dynamics. 
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Chapter 3 



Resummation in semi-inclusive DIS 



theoretical formalism 



Deep-inelastic lepton-hadron scattering (DIS) is one of the cornerstone processes to 
test PQCD. Traditionally, the experimental study of the fully inclusive DIS process 
e + A — > e + X, where A is usually a nucleon, and X is any final state, is used to 
measure the parton distribution functions (PDFs) for A. These functions describe the 
long-range dynamics of hadron interactions and are required by many PQCD calcula- 
tions. During the 1990's, significant attention has been also paid to various aspects of 
semi-inclusive deep inelastic scattering (SIDIS), for instance, the semi-inclusive pro- 
duction of hadrons and jets, e+A e+B+X and e+A e+jets+X . In particular, 
the HI and ZEUS collaborations at HERA, European Muon Collaboration at CERN, 
and the E665 experiment at Fermi National Accelerator Laboratory performed ex- 
tensive experimental studies of the charged particle multiplicity [^^ -|63fl and hadronic 
transverse energy fiows @,|6^ at large momentum transfer Q. It was found that 



some aspects of the data, e.g., the Feynman x distributions, can be successfully ex- 
plained in the framework of PQCD analysis |^, On the other hand, applicability 
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of PQCD to the description of other features of the process is limited. For example, 
the perturbative calculation in lowest orders fails to describe the pseudorapidity or 
transverse momentum distributions of the final hadrons. Under certain kinematical 
conditions the whole perturbative expansion as a series in the QCD coupling may fail 
due to the large logarithms discussed in Section O. 



To be more specific, consider semi-inclusive DIS production of hadrons of a type B. 
At large energies, one can neglect the masses of the participating particles. In semi- 
inclusive DIS at given energies of the beams, any event can be characterized by two 
energy scales: the virtuality of the exchanged vector boson Q and the scale qt intro- 
duced analogously to e~^e~ hadroproduction and Drell-Yan process (cf. Section [2^) . 
The scale qr is also related to the transverse momentum of B. The expansion in 
the series of as is justified if at least one of these scales is much larger than Aqcd- 
However, the above necessary condition does not guarantee fast convergence of per- 
turbative series in the presence of large logarithmic terms. If Aq^^ Q^, (It ^ Q'^ ■> 
the cross sections are dominated by the soft and collinear logarithms log'" {ItIQ"^) 5 
which can be resummed in the framework of the small-gr resummation formalism 



(Subsection |2.4.1| ). In the limit A^cd ^ It^ ^ It (photoproduction region) 
PQCD may fail due to the large terms log"* (Q^/o'j'), which should be resummed into 
the parton distribution function of the virtual photon ||6^. Finally, even in the region 
^QCD ^ <iT ^ encounter another type of large logarithms corresponding 

to events with large rapidity separation between the partons and/or the hadrons. 
This type of large logarithms can be resummed with the help of the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) formalism (Subsection |2.4.2| ). 

In this Chapter I discuss resummation of soft and collinear logarithms in SIDIS 
hadroproduction e + A — > e + B + X in the limit Ag^^ ^ Q^, ^ ■ The 
calculations are based on the works by Meng, Olness, and Soper [^,|3^, who analyzed 
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the resummation technique for a particular energy distribution function of the SIDIS 
process.^] This energy distribution function receives contributions from all possible 
final-state hadrons and does not depend on the specifics of fragmentation. 



Here the resummation is discussed in a more general context compared to p3|, p4|: 
namely, I also consider the final-state fragmentation of the partons. Using this for- 
malism, I discuss the impact of soft and collinear PQCD radiation on a wide class of 
physical observables including particle multiplicities. The calculations will be done 
in the next-to-leading order of PQCD. In the next Chapter, I compare the resumma- 
tion formalism with the HI data on the pseudorapidity distributions of the transverse 



energy fiow |^,^ and ZEUS data on multiplicity of charged particles ||6^ in the 
7*p center-of-mass frame. Another goal of this study is to find in which regions of 
kinematical parameters the CSS resummation formalism is sufficient to describe the 
existing data, and in which regions significant contributions from other hadroproduc- 



tion mechanisms, such as the BFKL radiation 56|, higher-order corrections including 



multijet production with |68| or without |69,70| resolved photon contributions, or 



photoproduction showering |71|, cannot be ignored 



3.1 Kinematical Variables 



I follow notations which are similar to the ones used in In this Section I 

summarize them. 

I consider the process 



e + A^e + B + X, 



(3.1) 



*The general features of the resummation formalism in semi- inclusive DIS were first discussed by 
J. Collins isl. 
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where e is an electron or positron, A is a proton (or other hadron in the initial state), 
-B is a hadron observed in the final state, and X represents any other particles in the 



final state in the sense of inclusive scattering (Fig. p. .21 ). I denote the momenta of A 
and B by and p^, and the momenta of the electron in the initial and final states 
by P and l'^. Also, is the momentum transfer to the hadron system, = l'^ — l'^. 
Throughout all discussion, I neglect particle masses. 

I assume that the initial electron and hadron interact only through a single photon 
exchange. Contributions due to the exchange of Z-bosons or higher-order electroweak 
radiative corrections will be neglected. Therefore, also has the meaning of the 4- 
momentum of the exchanged virtual photon 7*; q'^ is completely determined by the 
momenta of the initial- and final-state electrons. In many respects, DIS behaves as 
scattering of virtual photons on hadrons, so that the theoretical discussion of hadronic 
interactions can often be simplified by considering only the photon-proton system. 

3.1.1 Lorentz scalars 

For further discussion, I define five Lorentz scalars relevant to the process (|3.1|). The 



first is the center-of-mass energy of the initial hadron and electron \/SeA where 

SeA = {pA + If = 2pA ■ I. (3.2) 

I also use the conventional DIS variables x and which are defined from the mo- 
mentum transfer q'^ by 

g2 ^ _q2 ^ 21 . (3.3) 
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In principle, x and can be completely determined in an experimental event by 
measuring the momentum of the outgoing electron. 

Next I define a scalar z related to the momentum of the final hadron state B by 



Pb ■ Pa '^xpB ■ Pa 



The variable z plays an important role in the description of fragmentation in the final 
state. In particular, in the quark-parton model (or in the leading order perturbative 
calculation) it is equal to the fraction of the fragmenting parton's momentum carried 
away by the observed hadron. 

The next relativistic invariant gf. is the square of the component of the virtual 
photon's 4-momentum that is transverse to the 4-momenta of the initial and final 
hadrons: 

Qt = -QtQtt,, (3.6) 



where 



% =q^ - Pa Pb ■ 3.7 

Pa -Pb Pa- Pb 



As discussed in Subsection |2.4.1| , the momentum plays the crucial role in the re- 
summation of the soft and collinear logarithms. In particular, a fixed-order PQCD 
cross-section is divergent when — ^ 0, so that all-order resummation is needed to 
make the theory predictions finite in this limit. According to Eqs. ( ^75| , |377|) gf = if 
and only if = z {xp^ + g'^) . Hence the resummation is required when the final-state 
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hadron B approximately follows the direction of xpA + q- 

In the analysis of kinematics, I will use three reference frames. The most obvious 
frame is the laboratory frame, or the rest frame of the experimental detector. The 
observables in this frame are measured directly, but the theoretical analysis is com- 
plicated due to the varying momentum of the photon-proton system. Hence I will 
mostly use two other reference frames, the center-of-mass frame of the initial hadron 
and the virtual photon (hadronic cm., or hCM frame), and a special type of Breit 
frame which I will call, depending on whether the initial state is a hadron or a parton. 



the hadron or parton frame. As was shown in Ref. ||33|| , the resummed cross section 
can be derived naturally in the hadron frame. On the other hand, many experimental 
results are presented for observables in the hCM frame. These observables are not 
measured directly; rather they are reconstructed from directly measured observables 
in the laboratory frame. I will use subscripts h, cm and lab to denote kinemati- 
cal variables in the hadron, hCM or laboratory frame. Below I discuss kinematical 
variables in all three frames. 

3.1.2 Hadron frame 

Following Meng et al. ||33|,Q the hadron frame is defined by two conditions: (a) 



the energy component of the 4-momentum of the virtual photon is zero, and (b) the 
momentum of the outgoing hadron B lies in the xz plane. The directions of particle 



momenta in this frame are shown in Fig. 



In this frame the proton A moves in the +z direction, while the momentum transfer 
q is in the —z direction, and is 0: 

q^^ = (0,0,0,-Q), (3.8) 
p1,h = (3-9) 
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Figure 3.1: Geometry of the particle momenta in the hadron frame 



The momentum of the final-state hadron B is 



(3.10) 



The incoming and outgoing electron momenta in the hadron frame are defined in 



terms of variables ip and ip as follows 



Q 



^ — — (cosh'?/', sinh'?/' cos (yj, sinh-?/' sin (yj, —1) 



2 

Q 



i'l^ = — (cosh?/', sinh?/; cos sinh sin yj, +1) 



(3.11) 



Note that ip is the azimuthal angle of ih or i'l^ around the Oz-axis. ^ is a parameter 
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of a boost which relates the hadron frame to an electron Breit frame in which 



(g/2,0,0,-g/2). By (U) and (|3ll 



cosh'?/' = — — 1 = 1, 



y 



(3.12) 



where the conventional DIS variable y is defined as 



y 



xSgA 



(3.13) 



The allowed range of the variable y in deep-inelastic scattering is < ?/ < 1 (see 
Subsection p. 1.41 ); therefore >0. 

The transverse part of the virtual photon momentum has a simple form in the 
hadron frame; it can be shown that 



u / It c\ ^t-' 

<. = (-^,-gT,0,-- 



(3.14) 



In other words, is the magnitude of the transverse component of qt,h- The trans- 
verse momentum px of the final-state hadron B in this frame is simply related to gr, 

by 



(3.15) 



Also, the pseudorapidity of B in the hadron frame is 



r]h = -log(^tan 



dB,h 



1 (It 



(3.16) 



The resummed cross-section will be derived using the hadron frame. To transform 
the result to other frames, it is useful to express the basis vectors of the hadron frame 
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(T^, X^, Y^, Z^) in terms of the particle momenta [34|. For an arbitrary coordinate 
frame. 



X" 



Q ' 

e'^'^'^Z^TpX,, 



Q 



1 + ^ 

g2. 



(3.17) 



If these relations are evaluated in the hadron frame, the basis vectors T'^, X^^, F'' , Z^^ 
are (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1), respectively. 





a) Fragmentation of the quark current 



b) Fragmentation of tine target fiadron 



Figure 3.2: (a) In the current fragmentation region, the hadron-level cross section 
can be factorized into hard partonic cross sections aba, parton distribution functions 
Fa/A{^a, I^f)-, and fragmentation functions Ds/bi^b, I-^f)- (b) In the target fragmenta- 
tion region, the hadrons are produced through the mechanism of diffractive scattering 
that depends on "diffractive parton distributions" Ma^B/Ai^aXs, f^p)- 
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The limit of small Qt, which is the most relevant for our resummation calculation, 
corresponds to the region of large negative pseudorapidities in the hadron frame. 
Hence the resummation affects the rate of the production of the hadrons that follow 
closely the direction of the virtual photon. The region of negative rj^ is often called the 
current fragmentation region, since the final-state hadrons are produced due to the 
interaction of the virtual photon with the quark current. In the current fragmentation 
region, hadroproduction proceeds through independent scattering and subsequent 
fragmentation of partons. Therefore, in this region the hadron-level cross section 
gba can be factorized in the cross sections a^a for the electron-parton scattering 
e + a ^ e + 6 + X, the PDFs Fa/A(^a,M, and the FFs DB/b{^b,f^D) (cf. Figurelpa). 
The formal proof of the factorization in the current region of SIDIS can be found in 
3|,|7|. 



In the opposite direction rjh^ ^ {qt — ^ +00) contributions from the current frag- 
mentation vanish. Rather the produced hadron is likely to be a product of fragmenta- 
tion of the target proton, which moves in the +z-direction (cf. Eq. (p.9|) ). According 
to Eq. (|3.5|), such hadrons have z ~ 0. The target fragmentation hadroproduction is 
described by a different approach, which relies on factorization of the hadron-level 
cross section into cross sections of parton subprocesses and diffractive parton distri- 
butions Ma,B/A{^a,CB, fJ'p) (cf. Figure |3^b). These distributions can be interpreted 
as probabilities for the initial hadron A to fragment into the parton a, the hadron B, 
and anything else, and (s denote fractions of the momentum of A that are carried 
by the parton a and the hadron -B, respectively. The distributions Ma^B/Ai^aXs, fJ'p) 
(also called fracture functions) were introduced in Refs. ||73|,^ and used in [|67| ,|75H77| 
to describe various aspects of SIDIS with unpolarized and polarized beams. The fac- 
torization of cross sections in the target fragmentation region was formally proven in 
the scalar field theory and in full QCD |^^]. The recent experimental studies 
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of the diffractive scattering at HERA are reviewed in The detailed discussion 
of diffractive scattering and interesting models (8^ that are applied for its analysis is 
beyond the scope of this work. 



(I 



Y 



X 



Figure 3.3: Particle momenta in the hadronic center-of-mass (hCM) frame 



3.1.3 Photon- hadron center-of-mass frame 



The center-of-mass frame of the proton A and virtual photon 7* is defined by the 
condition pA.cm + qcm = 0. The relationship between particle momenta in this frame 



is illustrated in Fig. p.3| . As in the hadron frame, the momenta qcm and PA,cm in the 
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hCM frame are directed along the Oz axis. The coordinate transformation from the 
hadron frame into the hCM frame consists of (a) a boost in the direction of the virtual 
photon and (b) inversion of the direction of the Oz axis, which is needed to make the 
definition of the hCM frame consistent with the one adopted in HERA experimental 
publications. In the hCM frame the momentum of 7* is 

.L-(^^^,0,0,^^j, (3.18) 

where W is the hCM energy of the 7*p collisions, 



W^' = (PA + gf = gM^-l) >0- (3-19) 



Since all energy of the 7*p system is transformed into the energy of the final-state 
hadrons, W coincides with the invariant mass of the _B + X system. 
The momenta of the initial and final hadrons A and B are given by 

V%.^ = I 0' 0' I ' (3-20) 



P^Bcm^ [ EB,EBSmeB,cm:0,EBCOs9B,cm]: (3.21) 



where 



Eb = (3-22) 

COSeB,cm = , 2 - (3-23) 

The hadron and hCM frames are related by a boost along the ^-direction, so that 
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the expression for the transverse momentum of the final hadron B in the hCM frame 

is the same as the one in the hadron frame, 



Pt = ZQT. (3.24) 

Also, similar to the case of the hadron frame, the relationship between qt and the 
pseudorapidity of B in the hCM frame is simple, 

qj. = We-'^-"'. (3.25) 

Since the directions of the 2;-axis are opposite in the hadron frame and the hCM 
frame, large negative pseudorapidities in the hadron frame {qt 0) correspond to 
large positive pseudorapidities in the hCM frame. Hence multiple parton radiation 
effects should be looked for in SIDIS data at qr/Q ^ 1, or 



77,^ > In I I >2. (3.26) 



X 



The boost from the hadron to the hCM frame also preserves the angle ip between 
the planes of the hadronic and leptonic momenta, so that the momenta l^, l'^ of the 
electrons in the hCM frame are 

= I 4^ (^^' + V' + - Q') , I sinh cos 

-|sinh^sin<^, ^(^(iy2 + q2) + _ q2^ cosh^^ |; (3.27) 



C = I ^ [iW + Q') cosh ^-W' + Q'Y^ sinh ^ cos 
-^sinhV'sinv?, ^( -W^^-Q^ + (iy^-Q^)coshV ) )■ (3.28) 
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Finally I would like to mention two more variables, which are commonly used in 
the experimental analysis. The first variable is the flow of the transverse hadronic 
energy 



Et = Etotsm9cm, (3.29) 

where Etot is the total energy of the final-state hadrons registered in the direction of 
the polar angle 9cm- The measurement of Et does not require identification of indi- 
vidual final-state hadrons; hence Et is less sensitive to the final-state fragmentation. 
The second variable is Feynman x, defined as 

--%^-(l-|.V (3.30) 



In i ^M} Pb cm the longitudinal component of the momentum of the final-state 



hadron in some frame. For small values of qt, i.e., in the region with the highest 
rate, 

xf^z. (3.31) 
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Y 



X 




Figure 3.4: Particle momenta in the laboratory frame 
3.1.4 Laboratory frame 

In the laboratory frame, the electron and proton beams are collinear to the Oz axis. 
The definition of the HERA lab frame is that the proton (A) moves in the +z direction 
with energy Ea, and the incoming electron moves in the —z direction with energy E. 
The momenta of the incident particles are 



PAAab 



{Ea, 0, 0, E,, 



(3.32) 
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C=(E,0,0,-E). (3.33) 



We can use (|3 .21 , 13 .321) and ( p.33| ) to express the Mandelstam variable SeA in terms of 



the energies Ea, E in the lab frame: 

SeA = ^EaE. (3.34) 

The outgoing electron has energy E' and scattering angle 9 relative to the —z 
direction. I define the Ox-axis of the HERA frame in such a way that the outgoing 
electron is in the Ox2;-plane; that is, 

= {E', -E' sin e, 0, -E' cos 9) . (3.35) 

The four-momentum = — l'^ of the virtual photon that probes the structure of 
the hadron is correspondingly 

q^^^ = {E- E\ E' sin 9, 0, -E + E' cos 9) . (3.36) 

The scalars x and are completely determined by measuring the energy and the 
scattering angle of the outgoing electron: 

g2 = 2EE'(l-cos^), (3.37) 



EE'(l-cos^) 

" Ea[2E-E'{1 + cos9)\ ^^-^^^ 
Rather than working directly with E' and 9 (or and x), it is convenient to introduce 
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another pair of variables y and (3: 



g2 2E - E'il + cosO) 



and 



2xEa _ ^2EE'il-cose) 
^~ Q 2E-E'{l + cose)' ^ ' 



The variable y satisfies the constraints 



S, 



eA 



<y<h (3.41) 



where W is defined in the previous subsection. The relationship (|3.41|) can be derived 
easily by rewriting y as 

y = i- ^ =1 + ^, (3.42) 

^^eA ^eA 



where 



TeA = {pA - I'f. (3.43) 



Eq. ( p.41| ) follows from the geometrical constraints on TeA for the fixed invariant mass 



W"^ of the final-state hadrons: 

- SeA < TeA < 0. (3.44) 

The observed hadron (B) has energy Eb and scattering angle 9b with respect to 
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the +z direction, and azimuthal angle ips] thus its momentum is 



Pb lab = i^B, Eb sin 6b cos v^b, Eb sin 6b sin ipB, Eb cos 6b) 



(3.45) 



The scalars z and gf. depend on the momentum of the outgoing hadron: 



^Eb{1 - cos 6b) 

Q 



(3.46) 



2 '^EbEq 

Qt = 

z 



1 — cos 7 



(3.47) 



In Eq. ( |3.47| ) 7 is the angle between ps and xpA + q (cf. Fig.|3.4|); 



Q{l + {l-y)P^ 
2/3 



(3.48) 



is the energy component of xp'^ + q^. Define 6^ to be the polar angle of xp'^ + q- 



xP%iab + llab = Eo (1, sin 6,, 0, cos 6, 



(3.49) 



where 



cot^ = (3^1 -y. 



(3.50) 



The angle 7 in Eq. ( |3.47| ) can be easily expressed in terms of the angles 6^, 6b, and 

<fB, as 



C0S7 = cos^^^, cos + sin^* sin6'5 coslpb- 



(3.51) 
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Finally, the azimuthal angle of the lepton plane in the hadron frame (cf. 



Eqs. (|3.11| )) is related to the lab frame variables as 



cos if 



Q 



2qT \/l - 2/ 



1 , ?T 1 , 1 2 



(3.52) 



Figure [37^ shows contours of constant qr and in the plane of the angles 9b and (pB- 
The point = corresponds to = V^b = 0, in agreement with Eqs. (|3.47|J3.5l|) . 
According to these equations, qr depends ompB through cos ipB, which is a sign-even 
function of (pB- Thus each pair of gr, determines ipB up to the sign, so that the 
contours in Figure ^]5| are symmetric with respect to the replacement v^s — > — v^s. 
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^B,deg. 




6*6, deg 



x=0.01 

e, = 41° 

^B,deg. 




%,deg 



x=0.001 
6** = 156° 



Figure 3.5: The variables qt and (p as functions of the angles 6b, ^b- Solid lines 
are contours of constant qx for qr/Q ranging from 0.1 (the innermost contour) to 
3.0. Dashed lines are contours of constant for ip ranging from vr/lO to 37r/4. The 
contour = vr coincides with the 6'B-axis. The plots correspond to Ea = 820 GeV, 
E = 27 GeV, Q = 6 GeV, x = 0.01 (upper plot) and x = 0.001 (lower plot). 
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3.1.5 Parton kinematics 

The kinematical variables and momenta discussed so far are all hadron-level variables. 
Next, I relate these to parton variables. 

Let a denote the parton in A that participates in the hard scattering, with mo- 
mentum 

P^a = CaP% (3.53) 

Let b denote the parton of which is a fragment, with momentum 

(3.54) 

The momentum fractions and ^f, range from to 1. At the parton level, I introduce 
the Lorentz scalars x, z, qr analogous to the ones at the hadron level 

X — = — , 3.00 



^ Pb-Pa Z 

z^ = ^, (3.56) 

q-Pa 6 



Qt = (3-57) 



Here 5^ is the component of q'^ which is orthogonal to the parton 4-momenta and 
Pi 



qt-Pa^Qt-Pb^ 0. 
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Therefore, 



^ = <1'-P'a^-Pt^- (3.58) 

Pa -Pb Pa- Pb 

In the case of massless initial and final hadrons the hadronic and partonic vectors gf 
coincide, 

^ = (3.59) 



3.2 The structure of the SIDIS cross-section 

The knowledge of five Lorentz scalars SeA, Q, Qt, x, z and the lepton azimuthal 
angle (p in the hadron frame is sufficient to specify unambiguously the kinematics of 
the semi-inclusive scattering event e + A e + B + X . In the following, I will discuss 
the hadron cross-section das a, which is related to the parton cross-section ddi,a by 

daBA sr^ d^b . d^a ^ X dabai/J'F, IJ'd) 



j ^DB/b{^b:l^D) [ ^Fa/A{^a:l^F) 



dxdzdQ'^dq^dip J z ib ' Jx ia " dxdzdQ'^dq^dip 

(3.60) 

Here Fa/A{^a-i I^-f) denotes the distribution function (PDF) of the parton of a type a in 
the hadron A, and Ds/bi^b, I^d) is the fragmentation function (FF) for a parton type 
b and the final hadron B. The sum over the labels a, b includes contributions from all 
parton types, i.e., g,u,u,d,d, . . . . In the following, a sum over the indices i,j will 
include contributions from active flavors of quarks and antiquarks only, i.e., it will 
not include a gluonic contribution. The parameters /ip and /id are the factorization 
scales for the PDFs and FFs. To simplify the following discussion and calculations, I 
assume that the factorization scales iip-, I^d and the renormalization scale ii are the 
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same: 



f^F = f^D = fJ'- (3.61) 

The analysis of semi-inclusive DIS can be conveniently organized by separating 
the dependence of the parton and hadron cross-sections on the leptonic angle ip and 



the boost parameter ip from the other kinematical variables x, z, Q and qt |§3 



This separation does not depend on the details of the hadronic dynamics. Following 



3^ , I express the hadron (or parton) cross-section as a sum over products of func- 
tions of these lepton angles in the hadron frame Ap{ip,ip), and structure functions 
fVBAix,z,Q^,q^) (or fVbaix,z,Q^,q^,ij,), respectively): 



"^""^^ J2 'VsAix, z, Q^ q'r)A,{^, y^), (3.62) 



dxdzdQ'^dq^dip 



The coefficients ^Vba (or ''Ha) of the angular functions Ap{ip, ip) are independent of 
one another. 

At the energy of HERA, hadroproduction via parity- violating Z-boson exchanges 
can be neglected, and only four out of the nine angular functions listed in ||3^ con- 
tribute to the cross-sections (|3.62| - |3.63|) . They are 

Ai = 1 + cosh^ if)^ ^3 = — cos Lp sinh 2-?/', 

A2 = -2, = cos 2(p sinh^ ijj. (3.64) 
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a) 



b) 



c) 





d) 



e) 



f) 



Figure 3.6: Feynman diagrams for semi-inclusive DIS: (a) LO; (b-d) NLO virtual 
diagrams; (e-f) NLO real emission diagrams 

Out of the four structure functions, ^Ha for the angular function Ai = 1 + cosh^ 
has a special status, since only ^Vba receives contributions from the lowest order of 
PQCD (Figure |3^a). At 0{as), only the contribution to the ^Ha structure function 
diverges in the limit qx ^ 0. 

3.3 Leading-order cross section 

Consider first the 0{ag) process of the quark-photon scattering (Fig.|0|a). This 
process contributes to the total rate of SIDIS at the leading order (LO). There is no 
LO contribution from gluons. Due to the conservation of the 4-momentum in the 
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parton-level diagram, at this order 



Pb=Pa + (l ■ 



(3.65) 



This condition and Eqs. (|3.58|J3.59|) imply that 



Qt = -Qt ■ qt = 0. 



(3.66) 



Also the longitudinal variables are[J] 



X = Z =1, 



(3.67) 



so that the momentum of the final-state hadron B is 



(3.68) 



Since both quarks and electrons are spin-1/2 particles, the LO cross section is 
proportional to 1 + cos\^ if) = Ai^ifj^ip) (Callan-Gross relation [8^). Hence the LO 
cross section is 



da, 



dxdzdQ'^dq^dip 



LO "^e^ 



^KQt) 5ijeMl - x)5[l - z), 



(3.69) 



'''To obtain Eqs. (|3.67|) , con sider for insta nce, Eq. ( 3.65 ) in the Breit frame for v ~ Q and 3. By 
using explicit expressions (3. £- 3. 10 3.53 , 3.54 ) for the parton and hadron momenta at qr = 0, we find 



2 \x 



0, 



Eqs. (3.67) are solutions for this system of equations. 
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where qr = (0, qt-, 0, 0) in the hadron frame. In Eq. (p.69|) the parameter ctq collects 



various constant factors coming from the hadronic side of the matrix element, 



The factor Fi that comes from the leptonic side is defined by 

= Y^- (3-n) 

Cj denotes the fractional electric charge of a participating quark or antiquark of the 
flavor j; Cj = 2/3 for up quarks and —1/3 for down quarks. 

The LO cross section (|3.69|) does not explicitly depend on Q^, but rather on x 



and z. This phenomenon is completely analogous to the Bjorken scaling in completely 
inclusive DIS |||], i.e., independence of DIS structure functions from the photon's vir- 
tuality . Just as in the case of inclusive DIS, the scaling of the LO SIDIS cross 
section is approximate due to the dependence of PDFs and FFs on the factorization 
scale /i. This scale is naturally chosen of order Q, the only momentum scale in the LO 
kinematics. When /i ~ Q varies, the PDFs and FFs change according to Dokshitser- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) differential equations (|2.29| , |2.31|) . By solv- 
ing the DGLAP equations, one sums dominant contributions from the collinear ra- 
diation along the directions of the hadrons A and B through all orders of PQCD. 
Formally, the scale dependence of the PDFs and FFs is an 0(0:5) effect, so that it is on 
the same footing in the LO calculation as other neglected higher-order QCD correc- 
tions. By observing the dependence of the LO cross section on /i, we can qualitatively 
test the importance of such neglected corrections. One finds that this dependence is 
substantial, so that a calculation of 0{as) corrections is needed to reduce theoretical 
uncertainties. Let us now turn to this more elaborate calculation. 
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3.4 The higher-order radiative corrections 

The complete set of 0{as) corrections to the SIDIS cross section is shown in 
Figs, p .Gib- f. These corrections contribute to the total rate at the next-to-leading 
order (NLO). 

At this order, one has to account for the virtual corrections to the LO subprocess 
9 7* — > 9 (Figs.|0|b-d), as well as for the diagrams describing the real emission 
subprocesses Q 'j* Q g and g'y* — > qq, with the subsequent fragmentation of the 
final-state quark, antiquark or gluon (Figs. ^]6|e-f). The explicit expression for the 
0{as) cross section is given in Appendix 0. 

Due to the momentum conservation, the momentum of the unobserved final-state 
partons {e.g. the gluon in Fig. ^]6|e) can be expressed in terms of q^, p^, p'^: 



When there is no QCD radiation (p^ = 0), the momentum of b satisfies the leading- 
order relationships^ = Pa+q^, so that qx = 0. If gr/Q ^ 1, the perturbative parton- 
level cross section is dominated by the term with p = 1. In the limit q^/Q^ 0, but 
qr 7^ 0, ^Vba behaves as times a series in powers of and logarithms \n{q^/Q'^), 



where the coefficients t^^™ (x, z) are generalized functions of the variables x and 
Obviously, the coefficient of the order a| in Eq. (|3.73| ) coincides with the most diver- 
gent part of the C(a|) correction to the SIDIS cross section from the real emission 
subprocesses. This coefficient will be called the asymptotic part of the real emission 
correction to ^Vba at 0(a|). 



+ P^a-Pb- 



(3.72) 




(3.73) 
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Convergence of the series in (|3.73| ) deteriorates rapidly as qx/Q — ^ because of 
the growth of the terms (g^^) ln™(g|,/Q^). Ultimately the structure function ^Vha has 
a non-integrable singularity at = 0. Its asymptotic behavior is very different from 
that of the structure functions ^'^'"^Ha, which are less singular and, in fact, integrable 
at = 0. This singular behavior of ^Vha is generated by infrared singularities of 
the perturbative cross section that are located at qr = 0. Indeed, according to the 



discussion in Section|2^, the diagrams with the emission of massless particles generate 
singularities when the momentum of one of the particles is soft (pf 0) or 
collinear to the momentum of another participating particle {P1/3P2 = 0). The soft 
singularities in the real emission corrections cancel with the soft singularities in the 
virtual corrections. For instance, at 0{as) the soft singularities of the diagrams shown 
in Fig.|3.6|e-f cancel with the soft singularities of the diagrams shown in Fig.|3.6|b-d. 
The remaining collinear singularities are included in the PDFs and FFs, so that they 
should be subtracted from 

There exist two qualitatively different approaches for handling such singularities. 
The first approach deals with the singularities order by order in perturbation theory; 
the second approach identifies and sums the most singular terms in all orders of 
the perturbative expansion. In the next two Subsections, I discuss regularization of 
infrared singularities in each of these two approaches. 

3.4.1 Factorization of collinear singularities at 0{as) 

Let us begin by considering the first approach, in which singularities are regularized 
independently at each order of the series in as- The singularity in the 0{as) part of 
the asymptotic expansion (|3.73|) can be regularized by introducing a "separation scale" 



g| and considering the fixed-order cross section separately in the regions < qx < qx 



and qx > qx- The value of q^ should be small enough for the approximation ( p.73| ) 
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to be valid over the whole range Qt ^ Qt- 

The quantity gf plays the role of a phase space slicing parameter. In the region 



< Q't < Q't) can apply the modified minimal subtraction (MS) factorization 



scheme |g3| to take care of the singularities at qt = 0. In the MS scheme, the 
regularization is done through continuation of the parton-level cross section to n = 
4 — 2e, e > dimensions [^]. The n-dimensional expression for the 0{as) part of 



the asymptotic expansion (|3.73|) of ^Vba(x,^,Q'^,q^) is 



ba 



llql,0{as) ~ \ z 



4-n 



Q-qF; as 1 



5(1 - z) {P«(x) + P(,i)(x)} + {P«(z) + P(;)(z)} 5{l - x) 



+ 25{l-z)5{l-x) 



Cf log % - ^Cf 
# 2 



0{ — ,qrp 



\ TT 



(3.74) 



Here the color factor Cp = {N^ — l)/{2Nc) = 4/3, Ai'c = 3 is the number of quark 



colors in QCD. The functions Pj;j\^) entering the convolution integrals in ( |3.74|) are 



the unpolarized 0{as) splitting kernels ||47[|: 



i+e 



-(1-2^ + 2^^), 



c. 



i + (i-0^ 



The "+"-prescription in Pqq\C,) regularizes Pgq'iO at ,^ = 1; it is defined as 



(3.75) 

(3.76) 
(3.77) 



The scale parameter //„ in ( p.74| ) is introduced to restore the correct dimensionality 
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of the parton-level cross section dai,al {,dxd'zdQ'^dqJ^~'^dLp) for n 7^ 4. The soft and 
collinear singularities appear as terms proportional to and 1/e when n — 4. The 
soft singularity in the real emission corrections cancels with the soft singularity in the 
virtual corrections. At 0{as), the virtual corrections (Fig. [3^b-d) evaluate to 



da. 



ba 



dxdzdQ"^ dqj^dip 



virt,0{as) 27r V r(l - e) e 

ddba 



X 



dxdzdQ'^dq^dip 



(3.78) 



LO 



where the LO cross section is given in Eq. (p.69|) . 



The remaining collinear singularities are absorbed into the partonic PDFs and FFs. 
When the partonic PDFs and FFs are subtracted from the partonic cross section 
da, the remainder is finite and independent of the types of the external hadrons. 
We denote this finite remainder as {da)hard- The convolution of {da)hard with the 
hadronic PDFs and FFs yields a cross section for the external hadronic states A and 
B. The "hard" part depends on an arbitrary factorization scale /i through terms like 
Pj^l^ \n{jj,/K), where Pab\0 ^^e splitting functions, and K is some momentum scale 
in the process. The scales fj, and fj,n are related as 

The dependence on the factorization scale in the hard part is compensated, up to 
higher-order terms in as, by scale dependence of the long-distance hadronic functions. 

After the cancellation of soft singularities and factorization of collinear singulari- 
ties, one can calculate analytically the integral of ( ^Vba)hard over the region 
< qr ^ Qt- At 0{as) this integral is given by 
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S\2 



hard '2t[S£j\ 



(3.79) 



The LO and NLO structure functions are 



lyNLO 
ba,j 



5(1 - Z)6{1 - x)6bjSja, 



(3.80) 



5\2 



I (5(1 - ^)5(1 - x)5bjSja 



+ In- 



(1)/ 

bj 



+ 5(1 - z)6b,^;J'\x) + ci;'^'\mi - ^)s,a. 



(3.81) 



The coefficient functions 4^'°"*''^'' (0 that appear in ^V^alj^ given by 



^out(l) 



(0 = '^ijC'^ 



:(1 -0-25(1-0 



4f"(0 

{l)out 
^93 



2^(1 -0, 



(0 



(3.82) 

(3.83) 
(3.84) 



Now consider the kinematical region Qt > Qt^ where the approximation (|3.74| ) no 
longer holds. In this region, ( ^Vba)hard should be obtained from the exact NLO result. 
With this prescription, the integral over can be calculated as 



max (jfj. 



ddba 



dxd'zdQ'^dqj.dip 



Alii;, if)' 



{4? 



rpj p max Qj^ 

dqU'Vba) + / 

V /hard J(q^Y 



dqr ( ^Vba 



hard 



4 />max5|, 

+ VA,(V',<^) / dq^^l^Vba 

p=2 ^ 



(3.85) 



where maxg|, is the maximal value of q^ allowed by kinematics. The first integral 
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on the right-hand side is calculated analytically, using the approximation (|3.79| ); the 



second and third integrals are calculated numerically, using the complete perturbative 
result of the order 0{as). The numerical calculation is done with the help of a Monte 
Carlo integration package written in the style of the programs Legacy and ResBos 
used earlier for resummation in vector boson production at hadron-hadron colliders 



3.4.2 All-order resummation of large logarithmic terms 

A significant failure of the computational procedure in ( |3.85| ) is that it cannot be 
applied to the description of the gr-dependent differential cross sections. Indeed, 
the cancellation of the infrared singularities is achieved by integration of the cross 
section over the region < < Q'f- However the shape of the qt distribution is 
arbitrary and depends on the choice of the parameter gf that specifies the lowest qr 
bin < qT < Qt- The fundamental problem is that the terms in (|3.73|) with small 
powers of as do not reliably approximate the complete sum in the region qx Q. 

This problem justifies the second approach to the regularization of the singularities 
at = 0, in which large logarithms in ( |3.73| ) and virtual corrections at gr = are 
summed to all orders. A better approximation for ^V^a at g^/Q <^ 1 is provided by the 
Fourier transform of a 6-space function Wba{b,Q,x,^, fi), which sums the dominant 
terms in (p.73| ) and virtual corrections through all orders of as'- 

'VUx, % Q\ ql,fi)^ = ^ [ -^e^9^-%,„(6, Q, x, % /x). (3.86) 

Here 6 is a vector conjugate to gr, and h denotes the magnitude of h. Hence at 



70 



all values of Qt can be approximated by 



(3.87) 



where ^Yha is the difference between the 0[ots) expression for ^V^a (cf. Appendix 0) 
and 0{as) asymptotic part (p.74|) , taken at n = 4. This difference is finite in the 
limit gT — ^ 0. 

The complete hadron-level resummed cross section can be obtained by including 
the finite parton structure functions for p = 2,3,4 and convolving the parton-level 
structure functions with PDFs and FFs (cf. Eqs. (|3.60[ - p.63D ): 



da 



BA 



dxdzdQ^dqj^dip 



aoFiA.iij,^) r d% 
SeA 2 J (27r) 



-e"^^-'WBA{b,Q,x,z) + YBA. 



(3.88) 



In this equation, the hadron-level 6-dependent form-factor WBA{b,Q,x, z) is the sum 
of convolutions of parton-level form-factors Wba{b,Q,x,z) with the PDFs and FFs: 



Fa/A- 



(3.89) 



Yba denotes the complete finite piece, 



Yba = 'Yba + Yl '"^BAApi^j, y.), 

p=2 



(3.90) 



where 



^YbA = Y^ DB/b ®'Yba^ Fa/A- 



(3.91) 
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The explicit expression for Yba is presented in Appendix 

At small b and large Q (i.e., in the region where perturbative dynamics is ex- 
pected to dominate) the general structure of WBA{b, Q, x, z) can be found from first 
principles |[r7|,pT|: 

WBA{b, Q, x,z)=J2 e?(^B/6 ® C^;*)(^, &)(q^ ® Fa/A){x, 6)e-^^^(^'«). (3.92) 



According to the discussion in Section 2.^, the form-factor Wba is the all-order sum 



of the large logarithms, which remain after the cancellation of soft singularities and 
factorization of collinear singularities. The soft contributions are included in the 
Sudakov function SBA{b,Q)- At small 6, SBA{b,Q) does not depend on the types of 
the external hadrons and looks like 



SBA{b, Q)= ^ AiasiJi), Ci) In + B{as{Ti), Ci, C2) ] = 5^(&, Q), 



(3.93) 



with 

k 



A{as{ji),C,) = Y.MCM'^] , (3.94) 



k=l 

00 / V 



S(«5(/i),Ci,C2) = Y.B,{C,,C2)[^^^ . (3.95) 

Contributions from the collinear partons are included in the functions 0"'{x,b,^) 
and C°'^^{'z,b, fi). These functions can also be expanded in series of 0:5/7?, as 

q^ix, 6, /i) = C-^'^ix, Ci, C,,^,b) , (3.96) 

k=0 \ ^ / 
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q;^\z,b,^^) = J2cT^'\^,c„c„^^b)[^] . (3.97) 



A:=0 



According to Eg. (|3.93|) , the integration in S^{b,Q) is performed between two 



scales Ci/b and C2Q, where Ci and C2 are constants of order 1. These constants 
also appear in terms proportional to 5(1 — x) or 5(1 — z) in the C-functions. The 
complete factor W{b,Q) is approximately independent from Ci and C2. In addition, 
the C-functions depend on the factorization scale /i that separates singular collinear 
contributions included in the PDFs and FFs from the finite collinear contributions 
included in the C-functions. To suppress certain logarithms in 0{as) parts of the C*" 
and C°"* functions, it is convenient to choose 

Ci = 2e-^^ = 60, (3.98) 
C2 = 1, (3.99) 
fi = bo/b, (3.100) 

where •je = 0.577215... is the Euler constant. 

Using our NLO results, I find the following expressions for the coefficients Ak{Ci), 
Bk{Ci,C2) and the C-functions: 



^1 = Cf, (3.101) 
To the same order, the expressions for the C-functions are 



Bi = 2CFlog ^ ^ M - (3.102) 
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• LO: 



Qout{0) 
jk 

c 



z,fib) 

m(0) _ 
39 ~ 



c 



out(0) 
93 



6jk6{l - z), 
= 0; 



(3.103) 
(3.104) 
(3.105) 



NLO: 



CZ^'\x,f^b) = ^x(l-x)-P«(x)log(^), (3.107) 

C;f'\z,f^b) = %-n^;H^)log(^). (3.109) 

In these formulas, the indices j and k correspond to quarks and antiquarks, and 
g to gluons. In Appendix |^ I show that the expansion of the integral over b in 
Eq. ( p.88|) up to the order 0{as), with perturbative coefficients given in Eqs. ( |3.101 



p.l09| ), reproduces the small-gr limit of the fixed-order 0{as) cross section discussed 
in Subsection |3.4.1 . 



Due to the crossing relations between parton-level SIDIS, vector boson production, 
and e"'"e~ hadroproduction, the C*"-functions are essentially the same in SIDIS and 
the Drell-Yan process; and the C°"*-functions are essentially the same in SIDIS and 
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e~^e~ hadroproduction. At NLO the only difference stems from the fact that the 
momentum transfer is spacelike in DIS and timelike in the other two processes. 
Hence the virtual diagrams Figs. |3^b-d differ by vr^ for spacelike and timelike g^. 
Correspondingly, Cjl^^^ for SIDIS does not contain the term (7r^/3)5(l — x), which is 
present in the C*^''^^ -function for the Drell-Yan process. Similarly, C°f^^ for SIDIS 
does not contain the term (7r^/3)5(l — z), which is present in the Cj^*'^^'' -function for 
e'^e~ hadroproduction. 

Up to now, I was discussing the behavior of the resummed cross-section at short 
distances. The representation (|3.92|) should be modified at large values of the variable 
b to account for nonperturbative long-distance dynamics. The authors of Ref. |^ 
suggested the following ansatz for Wba which is valid at all values of b: 

WBA{b, Q, x,z)=Y, e^iDB/b ® C°,f){z, b,){qi ® Fa/A){x, b,)e-^^\ (3.110) 

3 

Here the variable 




(3.111) 



serves to reproduce the perturbative solution ( p.92| ) at 6 ^ bmax-, with bmax ~ 
0.5 GeV^\ and turn off the perturbative dynamics for b > bmax- Furthermore, the 
Sudakov factor is modified, being written as the sum of the perturbatively calculable 
part S^{b^, Q) given by Eq. ( 3.93| ), and a nonperturbative part, which is only partially 
constrained by the theory: 

SBA{b, Q, X, z) = S^'ib,, Q, X, z) + SZ{h, Q, x, z). (3.112) 
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An explicit solution for the function S^^{b,Q,x, z) has not been found yet. 
Nonetheless, the renormalization properties of the theory require that the Q depen- 
dence in the nonperturbative Sudakov term be separated from the dependence on the 
other kinematical variables, i.e., 

SsAib, Q, X, z) = g^BAib, X, z) + g%\i^, x, z) log (3.113) 

vo 

with Qo ~ 1 GeV. The theory does not predict the functional forms of (?^]^(6, a;, z) 
and g^^\{b,x, z), so these must be determined by fitting experimental data. In ad- 
dition, if S^^{b,Q,x, z) indeed describes long-distance dynamics, it should vanish 
or be much smaller than S^{b,Q,x, z) in the perturbative region b < bmax- In the 
analysis of the experimental results, we may find that the fit to the data prefers a 
parametrization of S'^^ib) that is not small in comparison to the perturbative part 
of Wba at 6 < bmax- Such observation will be an evidence in favor of important 
dynamics that is not included in the 6-space resummation formula with coefficients 
calculated at the given order of PQCD. Therefore, this work uses an interpretation 



of S that is broader than its original definition in ||2T|. S will parametrize not 



only large-6 physics, but additional contributions to Wba at all values of b that are 
not included in the perturbative part of Wba- In the following parts of the thesis I 
will test whether the data are consistent with the assumption that these additional 
contributions are small in comparison to the perturbative part of Wba when b < bmax- 
Before ending this section, I would like to comment on a subtle difference be- 
tween and C°"*. While the initial-state coefficient functions Cl^^^\x,Ci,C2,b, ^) 
given in Eqs. ( |3.106| ) and ( p. 107] ) depend on the factorization scale /i through a fac- 



tor In [fib/ bo], the final-state functions C^^^^^\'z, Ci, C2, b, jj) given in Eqs. (|3.108D and 



( p.lPgp depend instead on In [/i6/(6o?)]. The additional term oc Inz in the func- 
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tions C^"*''^'*(?, Ci, C2, 6, /i) becomes large and negative when ? — > 0, so that it can 
significantly infiuence the 0{as) contribution at small values ofz. As a result, the re- 
summed total rate tends to be lower than its fixed-order counterpart for z < 0.1. This 
issue is discussed in more detail in Section|4^. Similarly, the 0{as) part of the NLO 
structure function ^V,j^j^ in ( |3.81|) depends on /i through a logarithm In [yU^/(%|)^]. 

The appearance of the additional terms oc In ? in the functions C°^'^^^ and 
refiects the specifics of separation of the 0{as) "hard" cross section {da)hard from the 



collinear contributions to the FFs in the MS factorization scheme. The easiest way to 
see the specific origin of the In z' terms is to notice that the dependence on the param- 
eter jj^n in the n-dimensional expression ( p.74| ) for ^Vha{x,^,Q'^,q^) comes through a 
factor {2n fin / zY^"- , rather than through a more conventional (27r/i„)^~'^. In its turn. 



2^ appears in (27ryU„/?)^~'^, because the MS'-scheme prescribes to continue to ra — 2 di- 
mensions the transverse momentum of the outgoing parton, rather than the vector 
Qt = Vt/z relevant to the resummation calculation. It is this factor that generates 
the /i-dependent logarithmic terms In [fib/ [hoz)] in the functions C'^^^^\z^ Ci, C2, &, /i) 
and ^V^^j*^. The C*^*^^^-functions do not include In? because they are evaluated along 
the direction z = 1 in the phase space. In contrast, nothing forbids such a term in 
the functions C^"*^^'', in which 'z can be anything between z and 1. Moreover, the 



In? terms are needed to reproduce MS expressions for 0(0:5) coefficient functions in 
completely inclusive DIS ||85| by integration of dabal {dQ"^ dxdzdq^) over qr and z. 
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3.5 Hadronic multiplicities and energy flows 

Knowing the hadronic cross-section, it is possible to calculate the multiplicity of the 
process, which is defined as the ratio of this cross-section and the total inclusive DIS 
cross-section for the given leptonic cuts: 



1 da 

datot/dxdQ'^ dxdzdQ'^dq^dip 



Multiplicity ^ /.^.^2.^...o2w.2.,„ - (3-114) 



Both the cross-section and the multiplicity depend on the properties of the final-state 
fragmentation. The analysis can be simplified by considering energy flows which do 
not have such dependence. A traditional variable used in the experimental literature 
is a transverse energy flow {Et) in one of the coordinate frames, defined as 

{Et)^, = —J2 d^sEr^ L. (3.115) 

This definition involves an integration over the available phase space $b and a sum- 
mation over all possible species of the final hadrons B. Since the integration over $b 
includes integration over the longitudinal component of the momentum of -B, the de- 
pendence of {Et) on the fragmentation functions drops out due to the normalization 
condition 

J2 I z DB/b{z)dz = 1. (3.116) 

B 

Instead of (Et), I will analyze the fiow of the variable z. This fiow is defined 



as 86 



dx dQ"^ dq^dip ^ Jz^i„ dx dz dQ"^ dq^dip 
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I prefer to use rather than {Et) because {Ex) is not Lorentz invariant, which 
complicates its usage in the theoretical analysisfj. Since qx is related to the pseudora- 
pidity in the hCM frame via Eq. (|3.25| ), and the transverse energy of a nearly massless 



particle in this frame is given by 

Et ^ Pt = zqr, (3.118) 

the experimental information on dHz/idx dQ^ dq^) can be derived from the hCM 
frame pseudorapidity (r/cm) distributions of {Et) in bins of a; and Q^. If mass effects 
are neglected, we have 

dxdQ'^drjcmdifi dxdQ'^dq^dip 

By the factorization theorems of QCD, the hadron-level z-flow can be written 
as the convolution of a parton-level z-Qow with the PDFs, 



dxdQ'^dqj.dip 



ix "''^^^^"'^ dxdQ'^dq^dip ^ 



Similarly to the SIDIS cross section, the z-flow can be expanded in a sum over the 
leptonic angular functions Ap{ip,ip): 

-^^^ = pUi,QWr,,)AMM, (3.121) 

where the structure functions ^Vzaix, Q^, fi) for the 2;-flow are related to the struc- 
ture functions PVba{x,'z,Q'^,qT, fJ.) for the SIDIS cross section by 

■f-The 2:-fiow Ez is related to the energy distribution function E calculated in |Q as T,z — 
{2xEa/Q'^)'^- Here Ea is the energy of the initial hadron in the HERA lab frame. 
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'V,a{x, ql, lj) = y2 I '^ba{x, z, Q^ q^, fi). (3.122) 
, Jo 



The resummed z-flow is calculated as 



- / .n, ,3.1.3) 

where 

W.{b, Q,x) = J2 e?Cr* e-^^^'-'^'^) {C% ® F,/^) (x, 6., (3.124) 
i 

As in the case of the resummation of hadronic cross sections, only the structure 
function ^Vza for the angular function Ai = 1 + cosh^?/' has to be resummed. 
The functions Cj^J in ( gl24| ) are the same as in (|3.92| ). The coefficient C°"* is 



Cr- = l + --in' (3.125) 



The parameter 6^,, given by (|3TT1|) with = 0.5 GeV"\ is introduced in (|T24D 



to smoothly turn off the perturbative dynamics when b exceeds bmax- The term 
in ( |3.123| ) is the difference between the complete fixed-order expression at 0(0:5) for 
dJ2z/ {dxdQ^dq^d(p) and its most singular part in the limit qx 0; that is, 



Y '^^^ ( '^^^ 1 f 12 ) 

^ dxdQ'^dq^dip I dxdQ'^dq^dip 



asym 



The asymptotic part calculated to 0{as) is 
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dxdQ^dqj^dLf 



asym 



SeA 2q^ 2n 



It 



(3.127) 



Similar to (p.ll2| ), the z-Qow Sudakov factor 5*^ is a sum of perturbative and 
nonperturbative parts, 



S,{b, Q, x) = 5^(6,, g, x) + 5f ^(6, g, x) 



(3.128) 



The NLO perturbative Sudakov factor is given by the universal x-independent 
expression ( |3.93| ). As in the case of SIDIS multiplicities, the renormalization group 
invariance requires that the dependence of S'^^on In Q be separated from the depen- 
dence on other variables: 



(6, g, x) = g^'^ (b, x) + g^'^ (b, x) log 



(3.129) 



In principle, the 2;-flow Sudakov factor Sz{b,Q,x) is related to the Sudakov factors 
SBA{b,Q,x, z) of the contributing hadroproduction processes e + A^e + B + X 
through the relationship 



SBA{b,Q,x,z) I 



{DB/b(»C^f){z,b,,fi). (3.130) 



In practice, the efficient usage of this relationship to constrain the Sudakov factors 
is only possible if the fragmentation functions and the hadronic contents of the final 
state are accurately known. I do not use the relationship (|3.130|) in my calculations. 
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3.6 Relationship between the perturbative and re- 
summed cross-sections. Uncertainties of the cal- 
culation 

In the numerical calculations, some care is needed to treat the uncertainties in the 
definitions of the asymptotic and resummed cross sections, although formally these 
uncertainties are of order (^((as/Tr)^, g^^). 

3.6.1 Matching 



The generic structure of the resummed cross-section (|3.88|) , calculated up to the order 

0((a5/vr)^), is 



cr. 



resum ~ '^W + "* • (3.131) 



In (|3.131|) , the VT-piece receives all-order contributions from large logarithmic terms 



oo , 2k— 1 

iE(-) (3-132) 



m=0 



The F-piece is the difference of the fixed-order perturbative and asymptotic cross- 
sections: 

^ — ^pert ~ '^asym- (3.133) 

In the small-gT region, we expect cancellation up to terms of order 0{a^^^ /ti^^^) 
between the perturbative and asymptotic pieces in (|3.133|) , so that the H^-piece dom- 
inates the resummed cross-section ( |3.131|) . On the other hand, the expression for 
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the asymptotic piece coincides with the expansion of the W-piece up to the order 
0{ag /vr^), so that at large qr the resummed cross-section (|3.131| ) is formally equal 
to the perturbative cross-section up to corrections of order O^ag"^^ /n^~^^). 

In principle, due to the cancellation between the perturbative and asymptotic 
pieces at small qt, and between the resummed and asymptotic piece at large qt, the 
resummed formula (Tresum is at least as good an approximation of the physical cross- 
section as the perturbative cross-section apert of the same order. However, in the NLO 
calculation at qt ^ Q it is safer to use the fixed order cross-section instead of the 
resummed expression. At the NLO order of as, the difference between the H^-piece 
and the asymptotic piece at large qt may still be non-negligible in comparison to the 
perturbative piece. In particular, due to the fast rise of the PDFs at small x, the 
resummed and asymptotic pieces receive large contributions from the small-x region, 
while the perturbative piece does not (see the next Subsection for details). Therefore, 
the resummed cross-section aresum may differ significantly from the NLO cross-section 
apert- This difference does not mean that the resummed cross-section agrees with the 
data better than the fixed-order one. At qr > Q, the NLO cross-section is no longer 
dominated by the logarithms that are resummed in Eq. (|3.131| ). In other words, the 
resummed cross-section (|3.13lD does not include some terms in the NLO cross-section 
that become important at qr ~ Q- For this reason, at qr > Q the resummed cross- 
section may show unphysical behavior; for example, it can be significantly higher the 
NLO cross section or even oscillate if the W^-term changes rapidly near the boundary 
between the perturbative and nonperturbative regions. 

As the order of the perturbative calculation increases, the agreement between the 
resummed and the fixed-order perturbative cross-sections is expected to improve. In- 
deed, such improvement was shown in the case of vector boson production [^], where 
one observes a smoother transition from the resummed to the fixed-order perturbative 
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cross-section if the calculation is done at the next-to-next-to-leading order. Also, at 
the NNLO the switching occurs at larger values of the transverse momentum of the 
vector boson than in the case of the NLO. 

Since the fixed-order result is more reliable at > Q, the switching from the 
resummed to the fixed-order perturbative cross-section should occur at ~ Q. How- 
ever, there is no unique point at which this switching happens. Similarly, it is not 
possible to say beforehand which of the two cross sections agrees better with the data 
in the region ~ Q. In SIDIS at small x, the NLO 2;-fiow underestimates the data at 
Qt ^ Q, while the resummed z-flow is in better agreement. Therefore, it makes more 
sense to use the resummed z-fiow in this region, without switching to the fixed-order 
piece. On the other hand, in the charged particle production one has to switch to the 
NLO cross section at ~ <5 in order to reproduce the measured pT-distributions. 



3.6.2 Kinematical corrections at qr ^ Q 

In this Subsection I discuss the differences between the kinematics implemented in 
the definitions of the asymptotic and resummed cross-sections, and the kinematics of 
the perturbative piece at non-zero values of gr- 

Let us first discuss the NLO approximation to the hadronic cross section ( |3.60| ). 
The integrand of the NLO cross section contains the delta-function 



\x J\z 



1 



xz6 



(3.134) 



which comes from the parton-level cross-section ( |A.2| ). Depending on the values of 
X, z,Q'^,q^, the contour of the integration over C,a and determined by (p.l34| ) can 
have one of three shapes shown in Fig. [37^a,b,c. For qx the integration proceeds 
along the contour in Fig. |3.7| a, and the integral in ( |3.60|) can be written in either of 
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two alternative forms 

das A d^a 



dxdzdQ^dq^dip ^(^j^,^ ia-x 



f d^a 

/ £ ^ ^ MBA{^a,^b;X,Z,Q^,q^,<^) 



MBA{,ia,ih]X,z,Q'^ ,ql,if), (3.135) 



where 



DB/b{^b)Fa/A{L) Yl 'f''-^^^ ^' ^t)M^^ v)- (3.136) 

a,o p=l 

The lower bounds of the integrals are 

{ia)min = + ^' (3-137) 

(efe)mm = + (3-138) 

with 



= ^y/xz. (3.139) 

Alternatively, the cross-section can be written in a form symmetric with respect 
to X and z, 

'^"^ MsA+ t ^MsA, (3.140) 



dxdzdQ'^dq^dif Jx+w — x J^^^ C,b — z 

where the integrals are calculated along the branches RP and RQ in Fig. p77| a, re 
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spectively. As 0, 

(^a)mm ^ -^j (^fe)mm ^ 

and the contour PRQ approaches the contour of integration of the asymptotic cross- 
section shown in Fig. ^77| d. The horizontal (or vertical) branch contributes to the 
convolutions with splitting functions in Eq. (|3.74| ) arising from the initial (or final) 
state collinear singularities, while the soft singularities of Eq. (|3.74| ) are located at the 
point = X, ^b = z. 

On the other hand, as qt increases up to values around Q, the difference between 
the contours of integration of the perturbative and asymptotic cross-sections may 
become significant. First, as can be seen from ( p.l40| ), in the perturbative piece 



and are always higher than x + w ov z + w , while in the asymptotic piece they 
vary between x oi z and unity. At small x (or small z) the difference between the 
phase spaces of the perturbative and asymptotic pieces may become important due 
to the steep rise of the PDFs (or FFs) in this region. Indeed, for illustration consider 
a semi-inclusive DIS experiment at small x. Let qx/Q = 0.5, z = 0.5, and x = 10""^; 
then X + w = 3.2 ■ 10^'^ ^ x = 10~^. In combination with the fast rise of the PDFs 
at small x, this will enhance the difference between the perturbative and resummed 
cross-sections. 

Second, for x ov z near unity, it could happen that x + w > 1 or z + w > 1, which 
would lead to the disappearance of one or two branches of the integration of the 
perturbative piece (Fig. |37^b,c). In this situation the phase space for nearly collinear 
radiation along the direction of the initial or final parton is suppressed. Again, this 
may degrade the consistency between the perturbative and asymptotic piece, since 
the latter includes contributions from both branches of the collinear radiation. For- 
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tunately, the x — z asymmetry of the phase space in SIDIS is not important in the 
analysis of the existing data from HERA, since it covers the small-x region and is less 
sensitive to the contributions from the large z region, where the rate of the hadropro- 
duction is small. 

The numerical analysis below includes a correction that imitates the phase space 
contraction in the low-x region. This correction is incorporated by replacing x in 
Eqs. (|3.74 p.88| ) by the rescaled variable 



x=^^±^x. (3.141) 



This substitution simulates the phase space contraction of the perturbative piece. At 
small Qt, the rescaling reproduces the exact asymptotic and resummed pieces, but at 
larger qx it excludes the unphysical integration region of C,a ~ x. 
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Figure 3.7: The contours of the integration over ^q, for (a,b,c) the perturbative 
cross-section; (d) the asymptotic and resummed cross-sections 
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Chapter 4 



Resummation in semi-inclusive DIS 



numerical results 



Despite the abundance of experimental publications on SIDIS, none of them presents 
dependence of SIDIS observables on the variable qx- Hence gr-distributions, which 
are sensitive to the multiple parton radiation, have to be derived from the published 
data on less direct distributions. The gr-distributions for some of the HERA data 
were reconstructed for the first time in |^^]. In this Chapter, I concentrate on the 
analysis of the ^T-distributions for the 2;— flow ( ^3.1171) 

1 dE, 



atot dxdQ'^dq^ ' 

which can be derived from published pseudorapidity distributions for the transverse 
energy flow in the hCM frame |^,^. I will also discuss several observables, including 
the average value of q^, that were measured in the production of light charged hadrons 
6g,|6|. 



Reconstruction of the gr-dependence reveals an interesting trend in the data: 
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namely, the average qx (or average increases rapidly when either x or z decreases. 
This trend is illustrated in Figs.^]l| and Figure |4]^ shows the average gf. in the 



charged particle production for several bins of a; and z aX 2S < < 38 GeV^. The 
procedure of reconstruction of is described in detail in Section As can be 
seen from Figureji?^, in the ZEUS data sample ((x) = 1.94-10"'^) is several times 
higher than {q^) at the same value of {z) and larger values of (x) in the E665 data 
sample ((x) = 0.07 — 0.29). (gf.) increases even faster when {z) decreases and (x) is 
fixed. For instance, at (x) = 1.94 ■ 10~^ increases from 3 GeV^ at {z) = 0.775 to 
82 GeV^ at {z) = 0.075. 

A similar trend is apparently present in the behavior of the quantity a/ {Qt^z) / {^z), 
which was derived from the data for the distributions d {Et) /drjcm published in |^ . 
This quantity is shown in Figure ^]2| as a function of and x.[| At each value 



of Q^, a/ {Qt'^z) I i^z) becomes larger when x decreases. Also, a/ {Qt'^z) / {^z) is 
roughly constant along the lines of constant y = Q'^/xSeA {i-e., the lines parallel to 



the kinematical boundary y = 1). Larger values of a/ {Qt^z) / {^z) at smaller x are 
the evidence of "broader" distributions dJ^z/dgr- In the subsequent Sections, I discuss 
this phenomenon in the context of the g^-resummation formalism. 

In this Chapter I assume that the angle (p is not monitored in the experiment, 
so that it will be integrated out in the following discussion. Correspondingly, the 
numerical results for dasA/ {dxdzdQ'^dq^) will not depend on terms in Eqs. (|3.62 - 
|3.63|) proportional to the angular functions ^3 and A4, which integrate to zero. The 
dependence on the azimuthal angle ip is discussed in more detail in Chapter H. 



*The distributions \/ (g^S^) / (S]^\ were derived by converting distributions d{ET) /drjcm in 
dYjz/dq^ with the help of Eq. ( 3.1191 ) and then averaging (^q^T.^) and (Ez) over the experimen- 
tal bins of qr- In each bin of qr, central values of and were used. This procedure provides a 
reasonable estimate for {q'^Y,^) / (Sz) if the experimental gr-bins cover all available range of qr- 
Figure 4^ shows ^/(gl^EiyTT^iy for the "low-Q" data set of from which satisfies this criterion. 
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Figure 4.1: The average as a function of x and z in the charged particle production 
at = 28 — 38 GeV^. The data points are extracted from published distributions 
vs. xp |§0|,3 using the method described in Section|4^. 
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4.1 Energy flows 
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4.1.1 General remarks 

As was discussed at length in Subsection |3.4.2| , knowledge of the resummed SIDIS 
cross section can be used to predict the pseudorapidity spectrum of the transverse 
energy flow in the hadron Breit frame or the hCM frame. It is advantageous to study 
the energy flows, because they are less dependent on the specifics of the final-state 
fragmentation of the scattered partons into the observed hadrons. I therefore start the 
presentation of the numerical results with the comparison of the resummation formal- 
ism to the experimentally measured pseudorapidity distributions for the transverse 
energy fiow in the hCM frame. 



I consider the data on d {Et) /drjcm which has been published in ||6^,|65|. I consider 
seven bins of x and Q from g (10 < {Q^) < 50 GeV^ 3.7 ■ 10"^ < {x) < 4.9 ■ 10"=^) 
and two sets of bins of x and Q from ||65| ("low-Q^" set covering 13.1 < (Q^) < 
70.2 GeV^ 8 x 10"^ < (x) < 7 x 10"=^ and "high-Q^" set covering 175 < {Q^) < 
2200 GeV^ and 0.0043 < (x) < 0.11). 

The experimental distributions d (E-r) /drjcm at a fixed value of W"^ = Q^{l — x)/x 
can be converted into the distributions dJ^^/dqT using Eqs. ( p.25| , p.ll^ ): 



= We-"^^^ (4.1) 



and 

dH, 1 d{ET) 



(4.2) 



dxdQ'^dqx q^ dxdQ'^drjcm 
The "derived" data for dJ^^/ {dxdQ'^dq^) can be compared with the resummed z-fiow 
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( P.123| ), which is calculated as 



dxdQ'^dq^ SeA 
where 



"^^^ aoFi (1 + cosh^ ij) I ^e^^-%,(6, Q, a;) + F,, (4.3) 



WM Q,x) = J2 e?Cr* e-^^^'-'Q'^) iC% ® F,/^) (x, (4.4) 
i 

The Sudakov factor Sz in Eq. ( [4. 4]) is 

^,(6,Q,x) = S^(6*,Q,x) + ^f^(6,g,x), 
where the perturbative part is given by Eq. (p.93|), and a realistic parametriza- 



tion of the nonperturbative part S^^{b,Q,x) can be obtained by comparison with 
experimental data at low and intermediate values of Q, especially with the measured 
pseudorapidity distributions at Q ~ 3 — 20 GeV. At high Q, we expect the data to 
be dominated by the perturbatively calculable parton radiation and be less sensitive 
to the nonperturbative effects incorporated in S^^{b,Q,x). According to the renor- 
malization group invariance argument, S^^ includes a part that is proportional to 
InQ: 

S^^{b,Q,x)=g^'\b,x)+g^'\b,x)\og^, (4.5) 

vo 

where the parameter Qo ~ 1 GeV~^ prevents InQ/Qo from being negative in the 
region of validity of PQCD. In the following analysis, I use two parametrizations of 
S^^{b,Q,x), which I will call parametrizations 1 and 2. 



Parametrization 1 was proposed in our paper |^ with D. Stump and 
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C.-P. Yuan based on the analysis of the data in Ref. ||6^: 



5f^(6,Q,a;) = !^g^'\x) + i(^(2)(6,Q)|^^ + (6,g))} , (4.6) 

where g^'^\b,Q)\^Y g^'^\b,Q)\^^^_ are Q-dependent terms in the nonper- 
turbative Sudakov factors in Drell-Yan process and e+e~ hadroproduction. The 
parametrization of the function g^'^\b, Q) in Eq. ( [4.6|) is suggested by the cross- 
ing symmetry between SIDIS, the Drell-Yan process and e^e~ hadroproduction. 
Due to this symmetry, the functions g^'^\b, Q) in these processes may be related 



as 33 



9^"\^^ Q)\ SIDIS = I Q)\oy + 9^"\b, Q) L,,_) • (4.7) 

If the relationship (|4.71) is true, then the function g^'^\b,Q) in SIDIS is com- 
pletely known once parametrizations for the functions g^'^\b, Q) in the Drell-Yan 
and g^'^\b,Q) in e~^e~ hadroproduction processes are available. In practice, the 
only known parametrization of the nonperturbative Sudakov factor in the e+e^ 



hadroproduction was obtained in Ref. |19| by fitting the resummation formula 



to the data at Q = 27 GeV. Most of the (Et) data from HERA correspond to 
significantly smaller values of Q, where the usage of the parametrization (1^ is 
questionable. In addition, the known parametrizations of the nonperturbative 
Sudakov factors for the Drell-Yan P^, and e~^e~ hadroproduction 

processes correspond to slightly different scale choices: 

C^ = 6o, C2 = 1 (4.^ 
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and 



Ci = 60, C2 = e-''l\ (4.9) 

respectively. Therefore, the known functions (7*^^-' (6) and g^'^'^\^^^_{jS) are 
not 100% compatible and in principle should not be combined to obtain g^'^\b) 
for SIDIS. In the numerical calculation, I have used the functions g^'^\b)\^Y 
from [2^ and g^'^\b)\^^^_ from ll^, despite the fact that g^'^\b)\^y was fitted 



to Drell-Yan data using a different C2 value than in g^'^\b)\^^^_. Explicitly, the 
Q-dependent part g^'^\b,Q) in Eq. ( [4.6| ) is 



,(^.(,.Q).V(o.48.og(|-,+5.32C.,og(A).og(||)]. (4.10) 



In Eq. (CT), the constants are Ci = 2e~^^, C2 = e"^/\Qo = 1 GeV. The 



variable 6* is given by Eq. (|3TTTD , with b^ax = 0.5 GeV"^ 



The functional form of g^^\b,x) in terms of b and x was parametrized as 



g(^\b, x) = (-4.58 + (4.11) 



where the numerical coefficients were determined by fitting the experimental 
data. These data cover a limited region of x and {W < < 50 GeV^, 
3.7 ■ 10^^ < X < 4.9 • 10^'^ ), so that the parametrization 1 should not be 
used away from this region. Also, the dependence of S^^{b,Q,x) on Q can- 
not be determined reliably using exclusively the data from Ref. |^], since all 
pseudorapidity distributions in this publication are presented in a small range 
of Q ~ 2 — 6 GeV. This circumstance motivated us to model the Q-dependent 
terms in the parametrization 1 by using the crossing relationship (|4.71) instead 
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of trying to find these terms from the comparison with the data. 

Parametrization 2 overcomes several shortcomings of the parametrization 1. 
The parametrization 2 was proposed in ||36|| , where the analysis of Ref. 



was 

repeated using the latest and more comprehensive data on the transverse energy 
flow [^]. From our analysis, we found that the data from Refs. |§4|,^ are con- 
sistent with the following representative parametrization of the nonperturbative 
Sudakov factor: 



^f^(6,Q,x) =6^ ( 0.013^^^ + 0.19 In^ + C ) , (4.12) 

X Qo 



where the parameter Qo is fixed to be 2 GeV to prevent InQ/Qo from being 
negative in the region of validity of PQCD, and where we set C = for reasons 
explained later. 

The HI Collaboration presented pseudorapidity distributions of the transverse 
energy flow for up to 2200 GeV^. However, the data at such high is 
rather insensitive to the nonperturbative dynamics because of the poor resolu- 
tion of the HI detector in the region of large and rjcm- Thus the HI data at 
very high is not informative about the Q^-dependence of S^^{b, Q, x) either. 
Fortunately, the HI Collaboration presented distributions in two bins at inter- 
mediate values of (Q^), namely (Q^) = 59.4 QeV^ and (Q^) = 70.2 GeVl To- 
gether with the data from Refs. @, ^] at lower values of Q, these distributions 
provide the first direct tests of the Q^-dependence of S^^{b,Q,x). Therefore 
the parametrization 2 of S^^{b,Q,x) includes a numerical value for the coef- 
ficient of InQ/Qo, which yields reasonable agreement with all of the analyzed 
data. The resulting value for this coefficient differs noticeably from its model 
expression in the parametrization 1. However, we should not draw too strong a 
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conclusion from this difference, because it might be caused by ambiguities in the 
separation of dependence and x dependence in the existing data. To draw a 
strong conclusion about the crossing symmetry model, experimental pseudora- 
pidity distributions in a larger range of x at intermediate values of Q^, as well 
as improvements in the knowledge of the nonperturbative Sudakov factor in the 
e+e^-hadroproduction will be needed. 



4.1.2 Comparison with the data 

The numerical results below were obtained using the parameters of the HERA electron- 
proton collider. The energies of the proton and electron beams are taken to be 
equal to 820 and 27.5 GeV, respectively. All calculations were performed using 
CTEQ5M1 parton distribution functions |^] and the parametrization 2 of the non- 
perturbative Sudakov function S^^ (Eq. (|4.12|) ), unless stated otherwise. The theo- 



retical results in Figs. |4.3| - |45| were obtained using the kinematical correction to the 



asymptotic and resummed cross-sections at non-zero g^, which was discussed in Sub- 



section |3]6]^. The factorization and renormalization scales of the perturbative and 
asymptotic pieces are all set equal to ^ = Q. The resummed piece was calculated 
using Ci = &07 C2 = 1, fi = bo/b, where bo = 2e^^'^ . 

In Fig. [4.3| , I present the comparison of the existing data from [|6^ in one of the 
bins of X and {{x) = 0.0049, (Q^) = 32.6 GeV^) with the NLO perturbative and 



resummed z-flows. Figure demonstrates two important aspects of the NLO qr 
distribution (dashed curve): namely, the NLO z-flow exceeds the data at small qr and 
is below the data at qx > Q- In fact, the deficit of the NLO prediction in comparison 
with the data at medium and large qx {qr > 5 GeV) is present in the entire region 
of X and that was studied. 



As I discussed in Section |3]^, one can trust the resummed calculation only for 
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reasonably small values of qr/Q- For large values of gr, the fixed-order perturbative 
result is more reliable. This means that the NLO resummation formalism will not 
give an accurate description of the data for qr ^ Q due to the small magnitude of 
the NLO perturbative z-flow in this region. 



0.3 






0123456789 10 

qi. GeV 



Figure 4.3: Comparison of the NLO perturbative and resummed expressions for the 



2- flow distribution with the existing experimental data from HERA ||6J|. The data 
is for (x) = 0.0049, (Q^) = 32.6 GeV^. The resummed curve is calculated using the 



parametrization 1 of S^^. CTEQ4M PDFs [M were used. 



99 




qr, GeV 



qr, GeV 



qi, GeV 



Figure 4.4: Comparison of the resummed z-Qow (solid curve) in the current region of 
the hCM frame with the data in the \ow-Q^ bins from Refs. [^] (filled circles) and ||6^ 
(empty circles). For the bin with (Q^) = 33.2 GeV^ and (x) = 0.0047, the fixed-order 
0(0:5) contribution for the factorization scale ^ = Q is shown as the dashed curve. 
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Figure 4.5: Comparison of the resummed z-Qow in the current region of the hCM 
frame with the data in the high-Q^ bins from Ref. [^]. For the bin with (Q^) = 
617 GeV^ and (x) = 0.026, the 0(0:5) contribution for = Q is shown as a dashed 



curve. 
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The excess of the data over the NLO calculation at large qx (cf. Figs. ^73| - |4.5|) can 
be interpreted as a signature of other intensive hadroproduction mechanisms at hCM 
pseudorapidities rjcm < 2. A discussion of the cross-sections in this pseudorapidity 
region is beyond the scope of this thesis. There exist several possible explanations 
of the data in this region, for instance, the enhancement of the cross-section due to 
BFKL showering or resolved photon contributions |^,|7T|. It is clear, however, 
that better agreement between the data and the theory, in a wider range of ?7cm, will be 
achieved when next-to-next-to-leading order contributions, like the ones contributing 



to (2+1) jet production ||69|, are taken into account. 



On the other hand, Figs. ^73| - fl75| illustrate that the resummed 2;-flow is in better 
agreement with the data, over a wide range of qrlQ-, but also lies below the data 
if qrlQ significantly exceeds unity. The better consistency between the resummed 
z-flow and the data suggests that the resummed z-flow should be used up to values of 
^Zt/Q ~ 1 — 4, «. e., without switching to the fixed-order expression. This procedure 
was followed in the derivation of our numerical results. 

Let us discuss the features of the data presented in Figs.|4]^ and|4]^. First, the data 
in the low-Q^ bins is significantly influenced by nonperturbative effects and therefore 
is sensitive to the details of the parametrization of S^^{b, Q,x). This feature can be 
seen from the abundance of data points around the maximum of the gr-distribution, 
where the shape is mainly determined by S^^{b,Q,x). Also, the low-Q^ data from 
HERA is characterized by small values of x, between 10""^ and 10~^. For the theory 



to be consistent with the data from Ref. |6J| in this range of x, the nonperturbative 
Sudakov factor must increase rapidly as x — 0, at least as 1/y/x. Such x-dependence 
is implemented in the parametrization 1 of S^^. In our newer analysis, we found that 
growth of S^^{b, Q,x) as 1/x at small x is in better agreement with the more recent 
data from |6^ 
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Second, the data in the high-Q^ bins of Fig. ^]5| shows a behavior that is qualita- 



tively different from Fig. |4.4| . In the region covered by the experimental data points, 
the Qt distribution is a monotonically decreasing function of gr, which shows good 
agreement with the resummed 2;-flow over a significant rangeQ of qt- In the region 
qT < 10 GeV, i.e., where the maximum of the distribution is located and where 
nonperturbative effects are important, the experimental gy-bins are too large to pro- 
vide any information about the shape of dTj^/dqT. Thus, as mentioned earlier, the 
published high-Q^ z-Qow data from Ref. ||65| is not sensitive to the dynamics described 
by the nonperturbative Sudakov factor S^^{b,Q,x). 

A third comment is that most of the high-Q^ data points in Fig.|4^ correspond to 
(x) > 10~^. If the resolution of the HI measurements at high were better in the 
small-gT region, then the high-Q^ data would also reveal the behavior of S^^{b, Q, x) 
at large x. But, as mentioned before, the published data in the high-Q^ bins are not 
very sensitive to the shape of the z-flow at small qr. Therefore it is not possible to 
impose any constraints on S^^{b, Q,x) at large values of x, except that it should be 
positive, S^^{b,Q,x) > 0. For this reason we have chosen the x-dependent part of 
S^^{b,Q,x) in the parametrization 2 such that S^^{b,Q,x) grows approximately as 
1/x as X — i> and is positive for all x. For the same reason, we chose C = in the 
parametrization 2. Although the most general parametrization of S^^{b,Q,x) can 
have C 7^ 0, the current data cannot distinguish between the parametrization 2 with 
C = and C 7^ 0, as long as the value of C is not very large. 

Finally, Fig.|4.6| shows the results of our calculation presented as the hCM pseudo- 
rapidity distributions of the transverse energy flow (Et) . This quantity is obtained by 

the transformation ( |3.119|) . The small-gy region, where the resummation formalism 
'''I point out once again that both the Ojas) and resummed z-fiow lie below the data at very 



large qr, in all bins of x and Q in Figures i.4 and 4.5 
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is valid, corresponds to large pseudorapidities. In this region, the agreement between 
our calculation and the data is good. At smaller pseudorapidities (larger gr), one 
sees the above-mentioned excess of the data over the perturbative NLO calculation. 
In the {Et) vs. r/cm plot, this excess is magnified because of the factor in the 
transformation 



4.1.3 How trustworthy is the resummed z-flow at large g^? 

As noted earlier, the 0{as) fixed-order 2;-fiow is much larger than the data in the 
region qr/Q ^ I and smaller than the data in the region qr/Q ^ ^- In the small-gr 
region, the resummed z-How is, by its construction, more reliable than the fixed-order 
result. In the large-gr region, the resummed z-fiow, with the kinematical correc- 
tion ( ^.1411) included, is also in better agreement with the data than the fixed-order 



calculation. But theoretically, the resummed z-fiow at large qr/Q is not absolutely 
trustworthy, because it does not include those parts of the fixed-order 2;-fiow that are 
subleading in the limit qr ^ 0, but which might be important at large qx- If the 
NLO result were in a good agreement with the data at large qx, it would be justified 
to consider it a more reliable prediction in this region. But since the 0{as) contribu- 
tion is systematically smaller than the data, higher-order corrections are presumably 
necessary in order to describe the region qx ^ Q reliably. 

A systematic approach for improving the theoretical description of the large-gr 
region would require inclusion of the complete terms in both the fixed-order 

and resummed 2;-fiows. But because such a calculation is not available, it might be 
beneficial to use the resummed 2;-fiow as a better theoretical prediction both in the 
region qr/Q ^ 1, where application of the resummation formalism is fully justified, 
and for up to several units of Q, where the resummed z-fiow agrees with the data 
better than the fixed-order one. Then the use of the resummed g-r-distributions of 

104 



the z-Qow will provide more reliable predictions for other observables relevant to the 
SIDIS process. 

As an example, resummation can improve the reliability of the theoretical predic- 
tion for the azimuthal asymmetry of the z-Qow. The 6-space resummation formalism 
affects only the coefficient ^Vza of the angular function Ai{'ilj,ip). This coefficient 
is the one that dominates the yj-integrated ^-flow in the small-gy region, where the 
energy flow is the most intense. On the other hand, the main goal of the mea- 
surement of angular asymmetries is to study structure functions other than ^V^a, 
e.g., those corresponding to the angular functions A^{;ip,{p) = — cos sinh 2?/^ and 
Ai{^p,(p) = cos 2(y9 sinh^ t/). By using a better approximation for the coefficient ^Vza, 
it is possible to measure the coefficients ^'"^V^a more reliably. Conversely, by knowing 
that the all-order resummation effects are important in the region of small qt and by 
concentrating on the region where Qt is of the order Q or larger, one may find angular 
asymmetries that are well approximated in the lowest orders of PQCD. The impact 
of resummation on the angular asymmetries is discussed in more detail in Chapter^ 
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Figure 4.6: The hCM pseudorapidity distributions of the transverse energy flow in 
the current fragmentation region. The data are from CTEQ4M PDFs and the 
parametrization 1 of S^^ were used. 
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4.2 Normalized distributions of charged particle pro- 
duction 



Let us now turn to the discussion of particle multiplicities. Although the resummation 
formalism, as outlined in Chapter ^ can describe the cross section for any massless 
final-state particle (provided that the fragmentation functions for this particle are 
known), in this Section I will concentrate on distributions of the charged particle 
multiplicity, defined as 

1 da(A + e-^ + e + X) , , 

^ -. (4.13) 



atot dQ 

Here 9 is some kinematical variable, such as the variable in Eq. ( p.6| ), the transverse 
momentum pt of the final-state charged particle in the hCM frame, or the Feynman 
variable xp, 

,cm I -I (It 



Our calculation assumes that the charged particles registered in the detector are 
mostly charged pions, kaons and protons. Therefore the cross section for charged 
particle production can be calculated using (|3.60| ) with the replacement of the frag- 
mentation functions Ds/bi^b^ lA by 



Dh±/b{ib,^l)= DB/b{^b,fi)- (4.15) 



P,P 



The fragmentation functions DB/b{,ib,lj) are known reasonably well only for ^5 > 
0.05 — 0.1 |[87|-|89|. Thus, the formalism presented here is applicable to the production 
of charged particles with sufficiently large energies, i.e., for z > 0.05. 
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Certain experimental distributions are readily available from the literature ||59|,|60 



62|,^, such as da/dpx, da/dxp, as well as distributions for the average transverse 
momentum {p'^). However, the "experimental" qr distributions must currently be 
derived from pseudorapidity distributions by using Eq. ( p.25| ). Although the distri- 
butions da/dpT and {p"^) are quite sensitive to resummation effects, they cannot be 
interpreted as easily as the distributions da/dqT, primarily because the distributions 
da/dpT and {p\) mix resummation effects at small values of qx with perturbative 
contributions from the region qr/Q ^ 1- The most straightforward way to study the 
effects of multiple parton radiation would be to consider the g-r (or pseudorapidity) 
distributions that satisfy the additional requirement z > 0.05 — 0.1 and that are or- 
ganized in small bins of and x. Unfortunately, such distributions have not been 
published yet. Although Ref. |^ presents distributions da{p + e + e + X) / drjcm 

for some values of x and , these distributions are integrated over the full range of 
z. Therefore, they are sensitive to the uncertainties in fragmentation functions, mass 
effects,f| and contributions from diffractive scattering. 

Because the experimental qr distributions are unavailable, we have decided to 
undertake a simpler analysis than the one presented for the energy flow. Our goal 

here is to understand how the multiple parton radiation could affect various aspects 

■^Our calculation assumes that all participating particles, including the final-state hadrons, are 
massless. Because of this assumption, the production of final-state hadrons with z = is allowed. 
However, in realistic SIDIS experiments there is a non-zero minimal value of z determined by the 
finite mass of the observed hadron. It follows from the definition (3.5) of z and Eqs. (B.2C, 3.21) for 
the initial and final hadron momenta in the 7*p cm. frame, that 

where 

P%,cm — ^B.cm +PB,cm- 



According to Eq. (4.16), the mass of the final-state hadron should be included if z ~ ijib/W ^ 
^qcd/W . Hence, our massless calculation is not suited for the analysis of the distributions da{A-\- 
e + e + X)/dr]cm from ||5|], which are sensitive to such small values of z. 
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of charged particle production. For this purpose we focused our attention on data 
from the ZEUS Collaboration which presents the charged particle multiplicity 
in a phase-space region characterized by the mean values (W) = 120 GeV, (Q^) = 
28 GeV^, and the additional constraint z > 0.05. These values of (W) and (Q^) 
translate into an average value of a; = 1.94 x 10~^. A simple model for nonperturbative 
effects at small qr will demonstrate that resummation describes qualitative features 
of this set of experimental data better than the fixed-order calculation. 

In all of the cases presented, the strategy is to compare the resummed multiplicity 
to that from the next-to-leading order calculation. In the numerical analysis, the 
multiplicity was calculated using the CTEQ5M1 PDFs |Q and the FFs from |^. For 
the resummed multiplicity, the "canonical" combination Ci = bo, C2 = 1, /i = bo/b 
was used. The NLO cross section was calculated according to Eq. (|3.85| ), for the 
factorization scale fi = Q. As explained in detail in Section |3.4.1| , the integration 
of the NLO term over qx is done separately over the regions < qr < qx and 
Qt > Q'ri where q^ is a particular type of a phase space slicing parameter. The 
final results should not depend on the exact value of g| provided that it is chosen in 
the region where the 0{as) part of the next-to-leading-logarithmic expansion ( ^3.731) 
approximates well the exact NLO cross section. In practical calculations, q^ cannot 
be chosen to be too small, because the numerical calculation becomes unstable due 
to large cancellations between the integrals over the regions < qx < q^ and qx > 
q^. The NLO prediction for the integrated charged particle multiplicity (Tchgd/o'tot 
at (W) = 120 GeV, (Q^) = 28 GeV^ is practically independent of q^ in the region 
1 ~ Q't ~ 2.5 GeV (cf. Figure ^TTI) . The NLO distributions shown in the subsequent 
Figures were calculated for g| = 1.2 GeV, which lies within the range of stability of 
a^^^'^/atot- As in the case of the 2;-flow, the resummed charged particle multiplicity 
may suffer from matching ambiguities at gr/Q ~ 1- 
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Figure 4.7: The dependence of the 0{as) prediction for the total charged particle 
multiplicity on the value of the separation scale gf . The calculation is done for 
{W) = 120 GeV, = 28 GeV^. 



In Section ^T], we found that the resummed z-flow is in better agreement with the 
experimental distributions than the NLO z-flow, for the whole range Q't/Q ^ 2 — 4. 
That result suggests that it might be preferable to use the resummed z-flow in the 
whole range qx/Q ^ 2—4 as a better theoretical prediction, until the C(a|) prediction 
for the z-flow in the region qr/Q ^ 1 becomes available. In the case of the charged 
particle multiplicity, the resummed cross section, which is calculated according to the 
formula 



da 



BA 



dxdzdQ'^dq^ 



nS, 



eA 



d% 

(2^ 



e''^^-^WBA{h,x,z,Q)+YBA. 



(4.17) 
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Figure 4.8: The distributions (a) vs. Xp and (b) (q^) vs. Xp for the charged 
particle multiplicity at (W) = 120 GeV, (Q^) = 28 GeV^. The experimental points 
for the distribution (p^) vs. xp are from Fig. 3c of Ref. |^. The "experimental" points 
for the distribution vs. xp are derived using Eq. (|4.19| ). The solid and dashed 
curves correspond to the resummed and NLO (/i = Q) multiplicity, respectively. 



overestimates the experimentally measured rate for the production of charged parti- 
cles with pt > 2 GeV. This discrepancy indicates that the resummed cross section in 
the region qr/Q ^ 1 is too high, so that switching to the perturbative cross section 
in this region is in fact required. Therefore, we have chosen to use the resummed 
cross section for < 5 GeV and switch to the next-to-leading order cross section for 
gT > 5 GeV. 

As in the case of the z-flow, the shape of the qp distribution for the charged particle 
multiplicity at small values of qp depends strongly on the unknown nonperturbative 
Sudakov factor S^^{b,Q,x, z). For the purposes of this study, we introduced a pre- 
liminary representative parametrization of the nonperturbative Sudakov factor for 
the fixed values of x = 1.94 x 10^^ and = 28 GeV^, i.e., the values that coincide 



111 



with the average values of x and Q in ||60[. This z-dependent parametrization is 



{b, Q2 = 28 GeV^ x = 1.94 x 10'^, z) = (o.l8 + 0.8^^^^^ . (4.18) 

Since the ZEUS Collaboration did not publish pseudorapidity distributions for the 
charged particle multiplicity {1/ atot)da / drjcm in bins of varying z^ we had to deduce 
information about the ^-dependence of S^^ from the less direct distribution of 
vs. xf presented in Fig. 3c of [6^. This distribution, known as a "seagull" for its char- 



acteristic shape (Fig.[4.8|a), can be converted into the more illustrative distribution of 
vs. Xf (Fig.^^). Since the major portion of the registered events comes from 
the region q^/W"^ <^ 1, or Xi? ~ 2;, a first estimate of the experimental data points 
for the distribution of {q^) vs. Xf can be obtained by assuming that 



Wt) - ^ - (4.19) 



where {xf) denotes central values of xf in each bin in Fig.[1.8|a.F| We refer to the 
resulting values as "derived data". 

Note that the shapes of {p^) vs. xf and {q^) vs. xf are quite different. The 



transformation from Fig.|4.8|a to Fig.|4.8|b shows immediately that the wing-like shape 
of the distribution of vs. xf should be attributed to a purely kinematical effect, 
namely an extra factor l/z"^ which is absent in the distribution of {q^) vs. xf- Once 
this extra factor is removed, we see from Fig.|4.8|b that {q^) increases monotonically 



and rapidly as z approaches zero. In other words, the qx distribution broadens rapidly 
when z decreases. This behavior is approximately realized by the simple z-dependent 

nonperturbative Sudakov factor S'^^{b,Q,x, z) given in Eq. (|4.18| ). 

Hn principle, a more accurate experimental distribution (g^) vs. xp can be determined by its 
direct measurement. 
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Figure 4.9: The dependence of the charged particle multiplicity on the transverse 
momentum pt of the observed particles in the hCM frame. The data points are from 
60|| . The solid and dashed curves correspond to the resummed and NLO multiplici- 



ties, respectively. 

The parametrization of S'^^(6, Q, x, z) was chosen to maximize the agreement be- 



tween the resummed distribution of (gj.) vs. xp and the "derived data" (cf. Fig.[4.8|b). 



Figure |48|b shows that the resummed calculation is in better agreement with the data 



points than the NLO expression. We have found it difficult to reproduce the rapid 
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growth of as Xi? ^ in either approach. In the future, it will be interesting to 
see how a more precise theoretical study will be able to explain adequately this rapid 
growth of in the region xp —>■ 0, assuming that the actual experimental data for 
the {q^) vs. xp distribution resemble the "derived data" discussed above. 

The resummation also significantly affects the pt dependence of the charged par- 
ticle multiplicity. In Fig.^jD] we present the distribution {l/atot)d<^/dpT- We see 
that resummation effects must be included to describe the shape of this distribution 
at ^ 1 GeV. Furthermore, resummation also improves the agreement between 
the theory and the experiment in the whole range of Pt. Through Eq. (|]24D, the 
improved description of the qt distribution in the small-gr region translates into a 
better agreement with the pr distribution in the whole range of pt- Just as in the 
case of the z-flow, the fixed-order calculation gives a rate that is too small compared 
to the data, which implies that higher-order corrections are important. Until the 
complete corrections are available, the resummation formalism, which already 

accounts for the most important contributions in the region of the phase space with 
the highest rate {i.e., at small qr), serves as a better theoretical prediction in the 
whole range of pt- 



Finally, Fig. [4.10| shows the Xi^-distribution for the charged particle multiplicity 
{1/ atot)do- / dxp. We see that both the resummed and fixed-order distributions are 
in reasonable agreement with the data and with earlier published theoretical results 
for the 0{as) XiT'-distributions |^. For the fixed-order multiplicity, we present two 
additional curves corresponding to different choices of the factorization scale /x in 
( p.60| ); the lower and upper dotted curves correspond to /i = O.SQ and 2Q, respec- 
tively. Note that the scale dependence of the NLO multiplicity increases when z ^ 0. 
Also note that the resummed multiplicity is significantly lower than the data in the 
two lowest bins of xp {{xf) = 0.075 and 0.125), but consistent with the NLO mul- 
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tiplicity within the uncertainty due to the scale dependence. Such behavior of the 
resummed multiplicity results from the dependence of the 0{as) coefficient functions 
C^^^^^\'z, Ci, C2, yu) on the additional term In z which was given in Eqs. ( |3.108|) and 
( p.l09| ) and discussed at the end of Subsection |3.4.2| . This negative logarithm domi- 



nates the C"**^^) -functions at very small values of Similarly, the integral ( |3.79|) of 



the NLO cross section over the lowest bin < g|- < (gf)^ depends on In ^through 
the terms 



as given in ( |3.81|) . Numerically, this dependence is less pronounced than in the re- 



summed cross section. For z < 0.1, the growing scale dependence of the multiplicity 
in the 0{as) calculation indicates that unaccounted higher-order effects become im- 
portant and are needed to improve the theory predictions. For example, including 
the coefficient C^} in the resummed calculation will be necessary to improve 

the description of the charged particle multiplicity in the small-2; region. 
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Figure 4.10: The dependence of the charged particle multiplicity on the Feynman 
variable xp in the hCM frame. The solid curve corresponds to the resummed mul- 
tiplicity. The dashed, lower dotted and upper dotted curves correspond to the NLO 
multiplicity calculated for fi = Q, 0.5Q and 2Q, respectively. 
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4.3 Discussion and conclusions 



The results in this Chapter demonstrate that multiple radiation of soft and collinear 
partons influence a large class of observables and can be described with the help of 
the CSS resummation formalism P^,P3[. Multiple parton radiation affects hadropro- 



duction in the current region of deep-inelastic scattering, i.e., for large pseudorapidity 
of the final-state hadrons in the photon-proton cm. frame. 

Although the resummation formalism needs further development, in particular in 
the procedure for matching the resummed curve to the perturbative result in the 
transition region, it already improves the agreement between the theory and the data 
and provides interesting insights about qualitative features of SIDIS. The formalism 
describes well the behavior of the transverse energy flows measured at HERA ||64|,|65| 
in the region of large hCM pseudorapidity rjcm > 3.0. At smaller pseudorapidities 
the NLO rate falls below the existing data. Evidently, this is a signature of the 
importance of the NNLO corrections, which were not studied in this paper. The 
resummation formalism describes the pseudorapidity distributions of the transverse 
energy flow more accurately than the NLO calculation; this formalism also has good 
potential to improve the description of various distributions of particle multiplicity. 

The presented analysis shows that the experimentally measured qt distributions 
for the energy flow broaden rapidly as a; — > 0. This rapid broadening of the gr dis- 
tributions can be realized if the nonperturbative Sudakov factor in the resummed 
energy flow increases as 1/x. Similarly, the qx distribution for the charged particle 
multiplicity broadens rapidly when z 0, which is consistent with the nonperturba- 
tive Sudakov factor increasing as z'^-"^. The SIDIS nonperturbative Sudakov factors 
at small values of x and z are therefore qualitatively different from the known non- 
perturbative Sudakov factors for vector boson production and e~^e~ hadroproduction. 



117 



which do not depend on the longitudinal variables at all. The rapid growth of the 
nonperturbative Sudakov factor in SIDIS might indicate that the ep collider HERA 
tests the resummation formalism in a new dynamical regime, which was not yet stud- 
ied at colliders of other types. In particular, the CSS formula adopted here assumes 
the usual DGLAP physics for the evolution of initial and final state partons in 
which the radiation of unobserved collinear partons is kx-ovdeved. The broadening 
of the qT distributions may be a result of the increasing importance of fcr-unordered 
radiation in the limit a; — ^ 0. The growth of the nonperturbative Sudakov factor 5*^^ 
as X decreases may be caused by the increase of the intrinsic transverse momentum 
of the probed partons due to such radiation. 

There are several theoretical aspects of the resummation formalism that can be 
clarified when more experimental data are published. Perhaps the largest uncertainty 
in the predictions of the resummation formalism comes from the unknown nonpertur- 
bative contributions, which in the 6-space formalism are included in the nonpertur- 
bative Sudakov factor S^^{b). I have presented simple parametrizations of 5*^^(6) 
for the transverse energy fiow (cf. Eqs. ( |4.6|J4.12| )) and charged particle multiplicity 
(cf. Eq. (|4.18| )). These parametrizations were found by fitting the resummed energy 
fiow and charged particle multiplicity to the data from Refs. ||6^, |65|] and Ref. [|60 



respectively. Experimental measurements outside the range of those data will make 
it possible to further improve these parametrizations and, hence, the accuracy of the 
resummation formalism. 

The most straightforward way to study S^^{b) is by measuring the variation of 
the qt spectra of physical quantities due to variations of one kinematical variable, 
with other variables fixed or varying only in a small range. For the energy fiow, it 
would be beneficial to obtain more data at x > 10~^, where the predictions of the 
resummation formalism can be tested more reliably, without potential uncertainties 
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due to the small- a; physics. Another interesting question is the dependence of the 
nonperturbative Sudakov factor on the virtuality Q of the vector boson. This de- 
pendence can be tested by studying the qT spectra in a range of Q with sufficient 
experimental resolution in the current fragmentation region. Finally, to study effects 
of multiple parton radiation on semi-inclusive production of individual hadrons, it 
will be interesting to see the qT spectra for particle production multiplicities with 
the additional constraint z > 0.05 ~ 0.1, i.e., in the kinematical region where the 
parametrizations of the fragmentation functions are known reasonably well. 
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Chapter 5 

Azimuthal asymmetries of SIDIS 
observables 



In a recent publication [^T] the ZEUS Collaboration at DESY-HERA has presented 
data on asymmetries of charged particle (h^) production in the process e + p 
e + h"^ + X, with respect to the angle ip defined as the angle between the lepton 
scattering plane and the hadron production plane (of /i^ and the exchanged virtual 
photon). This angle is shown in Figure |H7^. The azimuthal asymmetries, (cosf/)) and 
(cos 2ip) , as functions of the minimal transverse momentum Pc of the observed charged 
hadron in the hadron-photon center-of-mass (hCM) frame, are defined as 

r f^^ dip cos ncp , , ,^2,1 — T- 
(cosn^)(p,) = /° Z (5.1) 



da 

dxdzdQ'^dpTdip 



with n = 1,2. In terms of the momenta of the initial proton p^, the final-state 
hadron p'^, and the exchanged photon g'^, the variables in (|5.1| ) are = ~q^q^, 
X = Q^/2(p^ ■ g), and z = {pa-Pb) I {pa- q)- J d^ denotes the integral over x, 2;, Q^,Pt 
within the region defined by 0.01 < a; < 0.1, 180 GeV^ < < 7220 GeV^ 0.2 < 
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z < 1, and Pt > Pc- Nonzero (cos 2^9) comes from interference of the helicity +1 and 
— 1 amplitudes of the transverse photon polarization; and nonzero {cos(p) comes from 
interference of transverse and longitudinal photon polarization. 

More than 20 years ago it was proposed to test QCD by comparing measured 
azimuthal asymmetries to the perturbative predictions ||9^. However, it was also 
realized that nonperturbative contributions and higher-twist effects may affect the 
comparison [[T^ , [QSHOSj . For example, intrinsic kx might be used to parametrize the 



nonperturbative effects ||9^, and indeed ZEUS did apply this idea to their analysis 
of the data. The relative importance of the nonperturbative effects is expected to 
decrease as px increases. Thus, the azimuthal asymmetries in semi-inclusive deep- 
inelastic scattering (SIDIS) events with large pr should be dominated by perturbative 
dynamics. 



By comparing their data to the PQCD calculation at the leading order in 05 [96- 
98| , the ZEUS Collaboration concluded that the magnitude of the measured asymme- 
try (cos ip) exceeds the theoretical prediction for Pc < 1 GeV, and (cos 2{p) is system- 
atically above the theoretical prediction for pc > 1.25 GeV. ZEUS also estimated the 
possible nonperturbative contribution, by introducing a transverse momentum kj^ of 
the initial-state parton in the proton, and similarly of the final-state hadron due to 
nonperturbative fragmentation. It was found that this nonperturbative contribution 
is negligible for {cos2{p). For (cosy?), the nonperturbative contribution can be sizable 
(up to 20%), but it is not large enough to account for the difference between the data 
and the 0(0:5) calculation at low Pc. Hence, it was suggested that the discrepancy in 
(cosf/)) may be caused by large higher-order corrections. 



From the comparison to the PQCD calculation at the leading order in as |96, 97 



the ZEUS Collaboration concluded that the data on the azimuthal asymmetries at 
large values of Pc, although not well described by the QCD predictions, do provide 
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clear evidence for a PQCD contribution to the azimuthal asymmetries. In this Chap- 
ter, the ZEUS data is discussed in a framework of QCD resummation formalism 
[TTI, [21|, |32|, |33|, |35|^7n that takes into account the effects of multiple soft parton emis- 



sion. The discussion targets two objectives. First, it is shown that the analysis of 
(cosy?) and (cos2(y9) based on fixed-order QCD is unsatisfactory because it ignores 
large logarithmic corrections due to soft parton emission. In addition, perturbative 
and nonperturbative contributions are mixed in the transverse momentum distri- 
butions, so that the presented data does not clarify the dynamical mechanism that 
generates the observed asymmetries. Second, I make two suggestions for improvement 
of the analysis of the ZEUS data. I show that perturbative and nonperturbative con- 
tributions can be separated more clearly in asymmetries depending on a variable qr 
related to the pseudorapidity of the final hadron in the hCM frame. I also suggest to 
measure the asymmetries of the transverse energy flow that are simpler and may be 
calculated reliably. I present numerical predictions for the asymmetries of transverse 
energy fiow. These predictions are the most important result in this Chapter. 

5.1 Large logarithmic corrections and resummation 

The resummation formalism applied here was discussed in Chapter |. It describes 
production of nearly massless hadrons in the current fragmentation region, where the 
production rate is the highest. In this region, transverse momentum distributions are 
affected by large logarithmic QCD corrections due to radiation of soft and collinear 
partons. The leading logarithmic contributions can be summed through all orders of 
PQCD [|^,|3^,|3^,|3^ by applying a method originally proposed for jet production in 



e'^e annihilation ||T^ and vector boson production at hadron-hadron colliders [21 



According to Eq. (|3.62|) , the spin-averaged cross section for SIDIS in a parity- 
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conserving channel, e.g., 7* exchange, can be decomposed into a sum of independent 
contributions from four basis functions Ap{ilj, if) of the leptonic angular parameters 



Here ifj is the angle of a hyperbolic rotation (a boost) in Minkowski space; it is 
related to the conventional DIS variable y, hj y = Q'^/xSeA = 2/(1 + cosh?/^). The 
angular basis functions are Ai = 1 + cosh^^/', A2 = —2, = — cos sinh 27/^, A4 = 
cos2y9 sinh^ t/). Of the four structure functions '^V, only and "^V contribute to the 
denominator of (p.l|), i.e., the (/j-integrated cross section. Of these two terms, ^V^ is 
more singular, and it dominates the rate. According to the discussion in Chapter^ 
the singular contributions in ^V" can be conveniently explored by introducing a scale 
qr related to the polar angle {6B,cm) of the direction of the final hadron (B) in the 



where ?7cm is the pseudorapidity of the charged hadron in the hCM frame. In the 
limit qT 0, the structure function is dominated by large logarithmic terms; it 



has the form qT'^YlT=ii^s/''^)^Ylm=o '^^''"^^^^"^ilT/Q'^)^ where ti^'^™) are some gener- 



alized functions. To obtain a stable theoretical prediction, these large terms must be 
resummed through all orders of PQCD. The other structure functions ^'^-'^I/ are finite 
at this order; they will be approximated by fixed-order 0{as) expressions. 

In Eq. ( ^.11) , the numerator of {cosip) or (cos 2(^9) depends only on the structure 
function or '^V, respectively. The measurement of (cosy^) or (cos 2^9) must be 




hCM frame: 



qr = Q\/l/x - 1 exp (-r^cm), 



(5.2) 
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combined with good knowledge of the yj-integrated cross section, i.e., the denominator 
of (|5.1|) , to provide experimental information on the structure function or ^V. 
Thus it is crucial to check whether the theory can reproduce the (/^-integrated cross 



section as a function of px before comparing the prediction for (|5.1|) to the data. But 
Figure ^]9| shows that the 0{as) fixed-order cross section is significantly lower than the 
data from in the range of Pt relevant to the ZEUS measurements. This difference 
signals the importance of higher-order corrections and undermines the validity of the 
0(0:5) result as a reliable approximation for the numerator of Eq. (|5.1|) . 

On the other hand, the resummation calculation with a proper choice of the non- 
perturbative function yields a much better agreement with the experimental data for 
the pT-distribution from |Q. One might try to improve the theoretical description 
of the ZEUS data using resummation for the denominator of Eq. ( |5.1| ). However, the 
resummation calculation for da / {dxdzdQ"^ dq^dip) in the phase space region relevant 
to the ZEUS data is currently not possible, largely because of the uncertainty in the 
parametrization of the nonperturbative contributions in this region. The impact pa- 



rameter (6-space) resummation formalism ^ includes a nonperturbative Sudakov 
factor , which contains the effects of the intrinsic transverse momentum of the initial- 
state parton and the nonperturbative fragmentation contributions to the transverse 
momentum of the final-state hadron (cf. Eq. (|3.112| )). Without first determining this 
nonperturbative factor, e.g., from other measurements, it is not possible to make a 
trustworthy theoretical prediction for the denominator of Eq. ( p.l|) and, hence, these 
azimuthal asymmetries. 
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Figure 5.1: Comparison of the 0{as) prediction for the ratio (cos y?) / (cos 2<^) with 
the ratio of experimentally measured values of (cosy)) and (cos2y9) from ||61|. The 
error bars are calculated by adding the statistical errors of (cosy^) and (cos 2^9) in 
quadrature. Systematic errors are not included. The theoretical curve is calculated 
for (x) = 0.022, (Q2) = 750 GeV^ using the CTEQ5M1 parton distribution functions 
90|| and fragmentation functions by S. Kretzer from [^]. 



The azimuthal asymmetries measured by ZEUS may also be sensitive to uncertain- 
ties in the fragmentation to in the final state. Indeed, the cross section in Eq. ([5T| ) 
includes convolutions of hard scattering cross sections with fragmentation functions 
(FFs), integrated over the range 0.2 < z < 1. Although the knowledge of FFs is 



steadily improving [p^-|89||, there is still some uncertainty about their z-dependence 
and flavor structure for the range of Q relevant to the ZEUS measurement. Therefore 
the most reliable tests of the theory would use observables that are not sensitive to 
the final-state fragmentation. The asymmetries {cos mp) would be insensitive to FFs 
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if the dependence on the partonic variable z were similar in the hard parts of the 



numerator and denominator of Eq. ( |5.1| ) , so that the dependence on the FFs would 
approximately cancel. It is shown in Appendix |^ that the partonic structure function 

, which dominates the denominator of (|5.1|), contains terms proportional to 
that increase rapidly as z decreases. However, the most singular terms in the partonic 
structure functions ^''^V are proportional to l/£ Therefore, the dependence on the 
FFs does not cancel in the azimuthal asymmetries. 

A curious fact appears to support the suggestion that the theoretical predictions 
for {cosTup) depend significantly on the fragmentation functions. While each of the 
measured asymmetries, {cosip) and {cos2ip), deviates from the 0{as) prediction, the 
data actually agree well with the 0{as) prediction for the ratio {cosip) / {cos2{p) , as 
shown in Fig.|5.1|. The error bars are the statistical errors on {cosip) and (cos 2(^9) 
combined in quadrature; this may overestimate the experimental uncertainty if the 
two errors are correlated. Since this ratio depends only on the numerators in Eq. (|5.1|) , 
which are less singular with respect to £ than the denominator, the dependence on 
the fragmentation functions may be nearly canceled in the ratio. The good agreement 
between the 0{as) prediction and the experimental data for this ratio supports our 
conjecture that the fragmentation dynamics has a significant impact on the individual 
asymmetries defined in Eq. (|5.1|) . 

The final remark about the azimuthal asymmetries in Eq. (|5.1|) is that the px (or 
Pc) distributions are not the best observables to separate the perturbative and non- 
perturbative effects. The region where multiple parton radiation effects are important 
is specified by the condition q^/Q"^ ^ 1- But the px distributions are smeared with 
respect to the qr distributions by an additional factor of z, because pr = z qx- Thus 
the whole observable range of pt is sensitive to the resummation effects in the region 
of qr of the order of several GeV. A better way to compare the data to the PQCD 
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prediction is to express the azimuthal asymmetries as a function of Qt, not px- Then 
the comparison should be made in the region where the multiple parton radiation is 
unimportant, i.e., for qr/Q ^ 1- 



5.2 Asymmetry of energy flow 

Next, I describe an alternative test of PQCD, which is less sensitive to the above 
theoretical uncertainties: measurement of the azimuthal asymmetries of the trans- 
verse energy flow. In the hCM frame, the transverse energy flow can be written as 

Unlike the charged particle multiplicity, the energy flow does not depend on the 
final-state fragmentation. According to the results in Chapter ||, a resummation cal- 
culation can provide a good description for the experimental data on the <y9-integrated 
ii^T-flow. A new class of azimuthal asymmetries may be deflned as 

f (i$ f^^ cos mp , ,i?Ji, dip 
{Et COS nifi) (qt) = 2. ^ • (5.4) 

J Jo dxdQ^dqj.d^^V 

The structure functions ^Ve^, for the i^^-flow can be derived from the structure func- 
tions for the SIDIS cross section using Eq. ( p.l22| ). Similar to the case of the par- 
ticle multiplicities, the asymmetries {Et cos (p) and {Et cos2{p) receive contributions 
from "^Vej. and ^Ve^, respectively. But, unlike the previous case, the denominator 
in Eq. ( |5.4|) is approximated well by the resummed Et-&ow. Thus these asymmetries 
can be calculated with greater confidence. 

Figure |5]2| shows the prediction for the azimuthal asymmetries {Ex cos ip) and 
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{Etcos2(^) as functions of Qt for (a) x = 0.0047, = 33.2 GeV^ in the left plots 
and (b) x = 0.026, = 617 GeV^ in the right plots. The asymmetries are shown 
in g^-bins that are obtained from the experimental pseudorapidity bins for the (p- 
integrated Et-How data from Ref. The upper x-axis shows values of the hCM 
pseudorapidity rjcm that correspond to the values of qx on the lower x-axis. For each 
of the distributions in Fig.|^, the structure functions '^Vej, and ^Vet were calculated 
at leading order in QCD, i.e., 0{as). The solid and dashed curves, which correspond 
to the resummed and 0{as) results respectively, differ because the structure function 
^Ve^ in the denominator of (|5.4|) differs for the two calculations. The resummed 
(/9-integrated Et-How is closer to the data than the fixed-order result, so that the 
predictions made by PQCD for the subleading structure functions ^Vej, and ^Ve^, 
will be confirmed if the experimental azimuthal asymmetries agree with the resummed 
distributions. 



The discussion in Section ^T] shows that in the region qt ^ Q the resummed 
(/^-integrated i^^-fiow is larger than the 0{as) prediction. This explains why the 
asymmetries for qt ^ Q are smaller for the resummed denominator than for the 
0{as) denominator. In the region Qt/Q ^ 1, the asymmetries are determined by 
the asymptotic behavior of the fixed-order and resummed partonic structure functions 
^Vej,. As Qt — > 0, the 0{as) structure functions ( ^VEj,)o{as)} ^^t) ^^^d "^Vej, behave 
as l/gf-, 1/gT and 1, respectively. Thus, asymptotically, the ratios ^'^Vej,/( ^VEj,)o{as) 
go to zero, although the qr distribution for the asymmetry {Ft cos ip) is quite large 
and negative for small, but non-vanishing qt (cf. Fig. ^.2|) . Resummation of ^Vej. 
changes the gy-dependence of the denominator, which becomes nonsingular in the 
limit qx 0. Consequently, the asymmetry {Et cos ip) with the resummed denomi- 
nator asymptotically grows as I/qt {i.e., in accordance with the asymptotic behavior 
of ^Vej,). Hence neither the fixed-order nor the resummed calculation for {Etcos(p) 
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is reliable in the low-g^ region, so that higher-order or additional nonperturbative 
contributions must be important at — 0. The asymptotic limit for the resummed 
{Et cos 2ip) remains finite, with the magnitude shown in Fig.|5]^. Since the magnitude 
of {Ej'Cos2ip) is predicted not to exceed a few percent, an experimental observation 
of a large asymmetry {Et cos 2ip) at small qx would signal the presence of some new 
hadronic dynamics, e.g., contributions from T-odd structure functions discussed in 
9|. 



Figure ^]2| shows that the predicted asymmetry {Et cos (p) {qt) at Qt/Q = ^ is 
about 1-2% for the resummed denominator, while it is about 2-4% for the 0{as) 
denominator. The asymmetry {Et cos 2ip) (qt) at Qt/Q = 1 is about 1.5-2% or 3.5- 
5%, respectively. Both asymmetries are positive for qt ~ Q. According to Fig.|5.2|a, 
the size of the experimental qt bins (converted from the t] bins in ||65|| ) for low or in- 
termediate values of is small enough to reveal the low-gr behavior of ^'"^Vej. with 
acceptable accuracy. However, for the high-Q^ events in Fig, the experimental 
resolution in qt may be insufficient for detailed studies in the low-g^ region. Nonethe- 
less, it will still be interesting to compare the experimental data to the predictions of 
PQCD in the region Qt/Q ~ 1, and to learn about the angular asymmetries at large 
values of and x. 

To conclude, the azimuthal asymmetry of the energy flow should be measured 
as a function of the scale qt. These measurements would test the predictions of the 
PQCD theory more reliably than the measurements of the asymmetries of the charged 
particle multiplicity. 
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Figure 5.2: Energy flow asymmetries {Ex cos ip) (qx) and {Ex cos 2{p) (qx) for (a) 
X = 0.0047, = 33.2 QeV^ and (b) x = 0.026, = 617 GeVl The Figure shows 
predictions from the resummed (solid) and the 0{as) (dashed) calculations. 
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Appendix A 



The perturbative cross-section, finite 
piece and z-flow distribution 



In this Appendix, I collect the formulas for the NLO parton level cross-sections 
ddba/ (dxd'zdQ'^dq^dip) , which were originally obtained in p8| . 

According to Eqs. ( p.60| , p.63D , the hadron level cross-section dasAl (dxdzdQ'^dq^dip) 
is related to the parton-level structure functions ^Vba as 



das A 



dxdzdQ'^dq^dif 



p=l 



(A.l) 



At non-zero qt, the parton cross-section receives the contribution from the real 
emission diagrams (Fig. ^]6^-f), so that ^V^a can be expressed as 



Vha - -7—^ 

AnSeA 71" 



Z 



''fba{x,z,Q ,qj,) 



(A.2) 



with the same notations as in Chapter 13. In this Equation, the qq structure func- 
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tions are 



''fjk{x,z,Q'^,qT) = 2CFXze]5jk "fjk, 

where 

'fit = 2(*/^i)=4i 
The ? 51 structure functions are 

where 

^fjg = 2 ( Vis) = 8; 
Finally, the g Q structure functions are 
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where 



J ( 



U ■ 



2(%)=4; 



(A.8) 



In Q), 



Qt = 



l-z' 



(A.9) 



The indices j and k correspond to a quark (antiquark) of a type j or /c, the index g 
corresponds to a gluon. 

The finite part Yba of the 0{as) cross section ( |A.1| ) is 



Yba 



cTpFi as 



p=l a,b 



(A.IO) 



For p = 1, the functions are 



^Rjk(x,z,Q^,qT) = S 



02 / V? 



fjkix, z, Q , qrp) 



+ 2Cf5{1 - z)5{l - x) ( log 2 



g2 2 



^Rjg{x,z,Q'^,qT) = S 



^-(l-i) fi-i 



fjg{x, z, Q , q'y) 



(A.ll) 



(A.12) 
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^Rgj{x,z,Q^,q^) = 5 



\x 



fgj{x, z,Q , qrp) 



„2 ^3 gq 



X) 



(A.13) 



For p = 2,3,4, 



,g2 



1 

- - 1 

X I \z 



Pha{x,z,Q\q^). (A.14) 



From (3.60), it is possible to derive the perturbative ^-flow distribution, 



E; zdz- 



das A 



dxdQ'^dqj.dif ^ Jzmin dxdzdQ'^dqj.dip 



a,b j p=l 



'ba{x,Z,Q ,qrp)Ap{llj,Lp). 



(A.15) 



It depends on the same functions ^ fbaix, Q^, g|.), with the parton variable z deter- 
mined by the 5- function in (|A.2|) , 



1 — a; 



{qllQ^-l)x + l' 



(A.16) 
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Appendix B 



0{ag) part of the resummed cross 
section 



In this Appendix, I demonstrate that the 0{as) part of the W-teiui in the resummed 
cross section ( [3.88| ) coincides with the small-gr approximation of the factorized 0{as) 
fixed-order cross section. Correspondingly the complete 0{as) part does not depend 
on the scales Ci/b and C2Q separating collinear-soft and collinear contributions to 
the resummed cross section. 



At qt <^ Q, the resummed cross section (|3.88D is dominated by the H^-term: 



das A 



dxdzdQ"^ dq^dip 



SeA 2 

apFi Ai{'il),ip) 

SeA 2Q2 



d% 

(27r 



,e''^"WBAib) 



Q 

■e—WBAr 



Q' 



fB.ll 



where 



(B.2) 
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According to Eq. (|3^ ), at 6 



(B.3) 



The perturbative Sudakov factor S'^and the C— functions in WBAib,Q,x, z) can be 
expanded up to 0{as) using Eqs. (|3.93| , |3.101hp^ 



dip' f C^O'^ 

A{as{jl). C) In + B{asm, C^] 



as (Ai , 2 



In^ + Si(a5(7i), Ci, C2) In ^^^^ ) + 



CI 



(B.4) 



iqi ® 6, /i) = Fa/A{x, /i) 1 + — 



05 



TC 



2 Ci 3 Ci 



F,/^)(a;,/x) - (In^Pj]) ® F,/^)(a;,/i)^ ; 
(DB/b ® C^f )(^, /i) = Dn/biz, fx) (^1 + ^Cf In^ 



(B.5) 



2 Ci 3 Ci 
H — in 



C^bn 2 C2&0 



where the functions cf^J^"' , cf^J"^^ are defined in Eqs. (|3.82|-|3.8^. In these equations the 



running of 05, which is an effect of is neglected. Inserting the 0{as) represen- 

tations ( p.4| - p.6| ) in Eq. (|B.3|) , we obtain the 0{as) expression for WBA{b,Q,x, z) : 



Wba{^,Q,x,z] 



0{as) j 



DB/j{z,ij)Fj/A{x,n) 



71 \2 bl 2 bl 



+ 



71 



zQ bo 



Pif){z,fx)^F,,A{x,f,) 
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Qbo 



(B.7) 



This expression does not depend on the constants Ci, C2, so that the only factorization 
scale in Eq. (|B.7|) is /i. The Fourier-Bessel transform of WBA^b, Q, x, z) to the g^-space 
can be realized by using relationships 



(2vr)^ 



(27r)2 



-iqx-h 



'^0 



5{qT 
1 



'^0 



2 

TT 



Ingf. 
9t 



(B.8) 
(B.9) 
(B.IO) 



where the "+"-prescription is defined as 



2tt p+ca 

d^^QT [/(gT)]+ giqr) = / dip grdgrfigx) ( giqr) - 9(0) ) • 

^0 



(B.ll) 



Hence the small-gr approximation for the (^(as) cross-section is 



da 



BA 



dxdzdQ^dqj^dif 



aoFiAi {ilj,ip) ^ 2 ^ 



S, 



eA 



X {5{qT)Fsix, z, Q, fx) + z, Q, qt, /i)) ,(B.12) 



where 



Fs{x, z, Q, /i) = DB/j{z, n)Fj/Aix, /i) + 



as 



pB/6®ci5)°"*)(^,/i)-pB/6®ln 



zQ 



+ Dn/j (z, /i) ( ® F^/a) (x, /i) - (In 



/i 



F/,')®F,/A)(x,/i)) )>; (B.13) 
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F+{x,z,Q,qT,fJ,) 



1 as 



2nQ^ TX 
2CFDB/j{z, n)Fj/A{x, n) 



Tt (It. 



0! 



X 



X {^{DB/b®P^f){z,^l)F,,A{x,^l) + DB/,{z.^l){p\^J ® (b.i4) 

The function Fs{x,z,fi) contributes at Qt = 0; it receives contributions from the 
leading order scattering, evolution of the PDFs between the scales /i and Q, evolution 
of the FFs between the scales fx and ^Q, and 0{as) coefficient functions c'^^^*", c*^i)°"*. 
The function F+{x, z,Q,qT, fi) is just a regularized asymptotic part ( p.74|) of the 
0{as) fixed-order structure function ^Vba (cf. Eq. ( p.74|) ). 



The 0{as) cross section (p.l2| ) can be integrated over the lowest bin of qx, 



< Q'l' < {qrY ^ Q^- The resulting integral is 



("t)^ _ da 



BA 



dxdzdQ'^dq'^dip 



S 



27r 



E 



Fs{x, z, Q, /x) + F[{x, z, ^) /^)^ • 



(B.15) 



where 



s 

F^{x,z, = 

— <t DB/j{z,fi)Fj/A{x,fl) 



Cp ( In^ 



3 In 



{4? 



+ In 



{4Y 



{DB/b®Plf){z,^i)F,,A{x,^i) 



+ DB/j {Z, li) {Pjl^ ® Fa/A){x, /i) 



(B.16) 



This expression agrees with Eqs. ( p.79| - |3.8l| ). Technically, this integration can be 
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easily realized by going back to the 6-space and using relationships 



Qj- S 

^ MqTb)qTdqT = ^JMrb): (B.17) 



r+oo V2 

/ Ji{ah)\u-^hdh = -ln(a2), (B.18) 

r+oo V2 

/ Ji{ab)\n^-^bdb = \n^{a^). (B.19) 
Jo % 
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